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LETTER OF THE EDITOR 


Liebe DERIVE- und Tl-Freunde, 

Ich wollte Ihnen gerne das vielseitige Programm 
der ACDCA Summer Academy vom 25.-28.8.1999 in 
Gosing, No vorstellen. Deshalb muss leider der Letter 
of the Editor sehr kurz ausfallen. 

Beachten Sie bitte das ausfuhrliche User Forum. Die 
Dateien konnen Sie wie letztlich von unserer Home 
page herunterladen. 

Mit den besten Gruflen und Wunschen fur einen 
schonen Sommer, Ihr 
Josef Bohm 


Dear DERIVE- and Tl-friends, 

I wanted to present the wide-ranging program of 
our ACDCA Summer Academy from 25 - 28 Au- 
gust. So you find only a very short '’letter” of the 
editor. 

Take note of the ample User Forum. You can down- 
load all files belonging to this issue from our home 
page. 

With my best wishes for a wonderful summer 
Josef Bohm 


ACDCA Summer Academy - Preliminary List of Talks 


Baumann, GER 
Biryukow/Fydorova, RUS 
Bohm, AUT 
Boieri, IT 

Ersoy/Adrahan, TUR 
Franzova, USA 
Graubner, GER 
Gremillion, USA 
Heugl, AUT 
Himmelbauer, AUT 
Keunecke, GER 
Knechtel, GER 
Kokolj-Volic, SLO 
Kutzler, AUT 
Lechner, AUT 
Leinbach, USA 
Lokar, SLO 
MaaB, AUT 
Magiera, POL 
Mahmudi, INDONESIA 
Mann, ISR 
Marlewski, POL 
Postel, GER 
Raggett, AUT 
Roanes, ESP 
Romanovskis, LAT 
Schuller, AUT 


Informationssicherheit durch kryptologische Verfahren 

Tzunami in DERIVE and TI-92 

Basic Skills versus Technology 

The LABCLASS Project in Italy 

Initiating a Project TI-92/DERIVE supported 

Using TI-92 in a Traditional Calculus and DE Course 

A Special Triangle 

TI Interactive! 

Important Elements of the Algebraic Competence 
TI-92 Programs in Maths Lessons 
Differential Equations 

Changes in Didactics/Methods using CAS-Calculators 
Exam Questions when using TI-92 or DERIVE 
Computers and Calculators as Pedagogical Instruments 
HIV and the Immune System - A Mathematical Model 
Using Computer Algebra to Extract Meaning from Parameters 
Exponential Growth 
Black Boxes and Teaching Mathematics 
Analytical Mechanics Problems with DERIVE 
Using MatLab/Maple in Solving Intervall Hull .... 

Are Circular Functions Trigonometric or Real? 

Lin. Discrete Least-Square Fitting Assisted by CAS 
A Tutorial System for High Schools 

Using the TI-92 and CBR in Physics and Applied Maths Teaching 
Truth-Tables in Propositional Multi-Valued Logics with DERIVE 
Kepler's Ellipses - From Compasses to CAS 
Curricular Interfaces 


Schwarz, GER Independent Classwork 

Shelby/Stoutemyer, USA DERIVE for Windows Version 5 - Rumor becomes Reality 

Sjostrand, SWE Computers and Mathematics at Lindalvs Gymnasium 

Sugeng, INDONESIA Maple and the Abstraction Process 

Tonisson, EST Step-By-Step Solution Possibilities in Different CAS 

Torres-Skoumal, AUT CAS in an International School Setting 

Weiskirch, GER Didaktische und Methodische Veranderungen im Leistungskurs 

Weller, GER Squaring the Circle and Leonardo 

Wumig, AUT Using the TI-92 in the 9 th Grade of Austrian Grammar Schools 
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Contributions: 

Please send all contributions to the Editor. 
Non-English speakers are encouraged to 
write their contributions in English to rein- 
force the international touch of the DNL. It 
must be said, though, that non-English 
articles will be warmly welcomed nonethe- 
less. Your contributions will be edited but 
not assessed. By submitting articles the 
author gives his consent for reprinting it in 
the DNL. The more contributions you will 
send, the more lively and richer in contents 
the DERIVE & TI-92 Newsletter will be. 


Next issue: September 1999 

Deadline 15 August 1999 


Preview: Contributions for the next issues 

3D-Geometry, Reichel, AUT 

Graphic Integration, Linear Programming, Various Projections a.o., Bohm, AUT 

A Utility file for complex dynamic systems, Lechner, AUT 

Examples for Statistics, Roeloffs, NL 

Quaternion Algebra, Sirota, RUS 

Various Training Programs for the Tl 

A critical comment on the "Delayed Assignation" :==, Kiimmel, GER 

Sand Dunes, River Meander and Elastica, The lighter Side Halprin, AUS 

Type checking, Finite continued fractions, Welke, GER 

Share Holder's Considerations using a CAS, Bohm, AUT 
Kaprekar's "Self numbers", GER 

Linear Programming - Graphic Solution on the Tl, Kirmse, GER 

Implicit Multivalue Bivariate Function 3D Plots, Biryukov, RUS 

LU-Decomposition, Morales, COL 

ODEs with Constant Coefficients, Fernandez, ARG 

Three Bodies Problem (DERIVE & Cyclone98), Fritsch, AUT 

and 

Setif, FRA; Vermeylen, BEL; Leinbach, USA; Speck, NZL; Biryukow, RUS 
Wiesenbauer, AUT; Aue, GER; Koller, AUT, Ibrahim & Cordoba, ESP 

Impressum: 

Medieninhaber: DERIVE User Group, A-3042 Wiirmla, D'Lust 1, AUSTRIA 
Richtung: Fachzeitschrift 
Herausgeber: Mag. Josef Bohm 
Herstellung: Selbstverlag 


The DERIVE-NEWSLETTER is the Bulle- 
tin of the DERIVE & TI-92 User Group. It 
is published at least four times a year with 
a contents of 44 pages minimum. The goals 
of the DNL are to enable the exchange of 
experiences made with DERIVE and the 
TI-92/89 as well as to create a group to 
discuss the possibilities of new methodical 
and didactical manners in teaching mathe- 
matics. 

As many of the DERIVE Users are also 
using the TI-92/89 the DNL tries to com- 
bine the applications of these modern tech- 
nologies. 

Editor: Mag. Josef Bohm 
A-3042 Wiirmla 
D’Lust 1 
Austria 

Phone/FAX: 43-(0)660 31 36 365 
e-mail: nojo.boehm@pgv.at 
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David Millward, USA 


davidmillward01ineone . net 


Hi, 

I am a registered user of Derive and bought all of the Derive texts but am having a problem solving 
difference equations using the RECUREQUN.MTH utility. (Now RecurrenceEquations.mth) 


In particular, the problem is with non-homogeneous equations such as 

u(n+l)-8u(n)+16u(n-l)=4n A 2 

and 

u (n) -6u (n-1 ) +9u (n-2) =3 A n 

Figures in brackets are subscript to denote the order of the term. 

CAN Derive be used to solve such problems ? If so, HOW !???? 

I do hope you can help me - it's driving me around the bend ! 

David Millward 


From Josef to Johann Wiesenbauer (my man for all cases!) 

As I am - unfortunately enough - not a specialist in difference equations I forward this request 
and hope that you can find an answer. Many thanks in advance, Josef 


Johann 's answer: 


To be honest, it doesn't exactly take "a specialist in difference equations" to solve the difference equa- 
tions mentioned above. I could serve as an example myself. Nevertheless, I had no problems whatso- 
ever when computing the following lines in DfW 4.1 1 (using RECUREQN.MTH for the very first 
time!) 

2 

2 2 ■ n 4- (3-n + 10- n + 11) 

LIN2_CCF (-§ , 16, 4-(n + 1) p n) = 2 -(cl + c2-n) + 

27 


2 ■ n 4- On + 10 -n + 11) 

ul(n) := 2 -(cl + c2-n) + 


27 


ul(n + 1) - S-Lil(n) + 16 ■ ul(n - 1) = 4-n 


n + 2 n 

LIN2_CCF (-6 , 9, 3 p n) = 3 ■ 


cl + c2-n + 


n 

2 j 


u2(n) := 3 


cl + c2-n + 


n n 


2 2 J 


u2(n) - 6 ■ u2(n - 1) + 9-u2(n - 2) = 3 


n 
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You see? It's simply a matter of evoking a library function. The online help for LIN2_CCF ( ) says 
“LIN2_CCF (p , q , r , x , cl, c2) simplifies to a general solution of the constant coefficient 
second order linear difference equation y(x+2) + /7y(x+ 1) + q-y(x) = r(x), in terms of two symbolic 
constants cl and c2.” 

i.e. just what the doctor prescribed! The only thing you have to take care of is the leading term of the 
difference equation which is u(n+ 2) and not u(n )\ ! Maybe it was this that drove David around the 
bend? Well, whatever... 

Cheers, Johann 

On my request there was also a very similar answer from David Stoutemyer, SWHH. Many 
thanks. 

See more 2015 comments on difference equations on page 57, Josef. 

Alfonso Poblacion, Valladolid, Spain alfonso6gauss .mat . eup . uva . es 

Dear Josef, 

Two comments about the last newsletter, #33. 

The following iterative algorithm 

"PIPOB.MTH" 

Y (n) : =ITERATES ( (1- (l-k_ A 4) A (1/4) ) / (1+ (l-k_ A 4) A (1/4) ) , k_, SQRT (2) - 1 , n) 

A (n) : =IF (n = 1, (1+ (Y (n) ) j (n+1) ) A 4 * ( 6-4 • #2 ) -2 A (2n+l) * (Y (n) ) j (n+1) * 
(l+(Y(n) ) | (n+1 ) + ( Y (n) ) ; (n+1) A 2) , 

(1+ (Y (n) ) | (n+1) ) A 4 *A (n-1 ) -2 A (2n+l) * (Y (n) ) ; (n+1) * 

( 1+ ( Y (n) ) | (n+1 ) + ( Y (n) ) ; (n+1) A 2) ) 

PIAPPROX (n_) :=1/A(n_) 

gives approximations of PI: approximating PIAPPROX(l), you have 8 exact decimals; 
PIAPPROX(2), 41; PIAPPROX(3), 171;PIAPPROX(4), 694; etc, and only in a few seconds (I used 
DFW). As I don't work with TI92, 1 don't know if these good results are also obtained with it, so any 
interested user could check it and tell us. I think that it should be faster than the seven hours of the 

Spigot algorithm. With this iterative method, you don't need to store previous digits either. (I include 

this algorithm in a chapter of a book for secondary levels that it is in press by now). 

There are other algorithms based on modular equations. 

On the other hand, you comment some implementations of the Koch curve in page 34. There is an- 
other one, extensive to any Iterated Function System: the one I explained in Plymouth, do you remem- 
ber? It can be seen in the last chapter (’’Construction of Fractal Sets”) of the book ’’Mathematical Ac- 
tivities with DERIVE”, edited in 1997 by Chartwell-Bratt. It also contains files to represent the fractal 
structure of some types of infinite series. 


Cheers, Alfonso 
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Preci sionDigi ts := 8 
PIAPPROX(l) 

3.1415926 

Preci sionDigi ts := 41 
PIAPPR0X(2) 

3.1415926535897932384626433832795028841971 

Preci sionDigi ts := 171 

PIAPPR0X(3) 

3.141592653589793238462643383279502884197169399375105820974944592307816406286208 

99862803482534211706798214808651328230664709384460955058223172535940812848111745 

028410270193 

Preci sionDigi ts := 694 
PIAPPR0X(4) 

3.141592653589793238462643383279502884197169399375105820974944592307816406286208 

99862803482534211706798214808651328230664709384460955058223172535940812848111745 

02841027019385211055596446229489549303819644288109756659334461284756482337867831 

65271201909145648566923460348610454326648213393607260249141273724587006606315588 

17488152092096282925409171536436789259036001133053054882046652138414695194151160 

94330572703657595919530921861173819326117931051185480744623799627495673518857527 

24891227938183011949129833673362440656643086021394946395224737190702179860943702 

77053921717629317675238467481846766940513200056812714526356082778577134275778960 

9173637178721468440901224953430146549585371050792279689 

“DERIVE ’s built in tt” 

7T 

3.141592653589793238462643383279502884197169399375105820974944592307816406286208 

99862803482534211706798214808651328230664709384460955058223172535940812848111745 

02841027019385211055596446229489549303819644288109756659334461284756482337867831 

65271201909145648566923460348610454326648213393607260249141273724587006606315588 

17488152092096282925409171536436789259036001133053054882046652138414695194151160 

94330572703657595919530921861173819326117931051185480744623799627495673518857527 

24891227938183011949129833673362440656643086021394946395224737190702179860943702 

77053921717629317675238467481846766940513200056812714526356082778577134275778960 

9173637178721468440901224953430146549585371050792279689 


(Compare with Peter Witthinrich's results in DNL#33, Josef) 


G P Speck, Wanganui, New Zealand 

You might find the following comments, in whole or in part, suitable for inclusion in the Forum sec- 
tion of the DERIVE Newsletter. 

I found the questions and discussion from Hugo, Jakobcic, Wiesenbauer and Stegenga concerning 
piecewise defined functions on pages 5 & 6 of DNL # 32 in the DERIVE User Forum quite interesting. 
All of the issues and shortcomings raised in this discussion can be resolved easily using the character- 
istic functions defined in my article Concentric Curve Shading in DNL# 3 1 . 

(The functions can be found in the file SHAD_UT.MTH on the Diskette of the Year 1998. Or 
you can download the file from SWHH's home page - MTH31. Look for the DUG-section. 
Josef) 

Now it seems that there are three issues raised concerning piecewise defined functions expressed via 
characteristic functions in the DNL#32 : 
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1) Johann Wiesenbauer rightfully points out that CHI ( a , x , b ) as defined is not a function since 
SIMPLIFYING, e.g., CHI (0,0,1) or CHI (0,1,1) gives ± 0.5 + 0.5. This deficiency is cor- 
rected in our CHOPEN ( x , a , b ) characteristic function; it has value 1 at each x in the open interval 
(a,b) and has value 0 elsewhere. It is essentially the real-valued counterpart of CHI ( a , x , b ) with 
restrictions, a < b, a ^ oo, and b ^ go, which pose no major problems with characteristic functions 
such as CHLESS (x , a) and CHGRT (x , a) defined. 

2) Johann Wiesenbauer suggests defining a function, say CHI2(a,x,b,c,d) which is 1 on the open 
interval (a,b), has value c at x = a, has value d at x = b, and is 0 elsewhere. This is, of course, trivial 
using our characteristic functions: 

CHI2 (a, x,b, c, d) :=CHOPEN (x, a,b) +c*CHPOINT (x, a) +d*CHPOINT (x,b) 

3) To display Allan Hugo's piecewise defined function using characteristic functions, first con- 
sider: 

G (x) : =x*CHOPENRT (x, 0, 1 ) +x A 2 *CHOPENRT (x, 1,3) + (2*x+3) *CHGRTEQ (x, 3) 
[G(-1),G(0),G(1),G(3)] = [0,0,1,9] 

G(x) is a characteristic function representation of Allan Hugo's function except that G is defined as 
0 for x < 0 while Hugo's function is undefined for x < 0. Johann, again rightfully, points out that 
this is not satisfactory and feels that resort to the use of one IF() seems to be necessary. However, 
even here, we can escape using an IF(). Consider: 

F (x) : =FLOOR ( SQRT (x-10) ) *CHLESS (x, 0) +G (x) 

Note that the FLOOR() function is undefined at any non-real complex number, so that 
FLOOR(SQRT(x-10)) is undefined at any x < 0. The following values of F illustrate that it is unde- 
fined at x < 0: 

[F (-100) ,F(-1) ,F(0) ,F(1) ,F(3) ] = 

[FLOOR (SQRT (110) *1) , FLOOR (SQRT (ll)*i),0,l,9] 

This does require, of course, that when we get FLOOR(c), where c is non-real complex, returned 
when we SIMPLIFY F(r), we interpret this as meaning that F is undefined at x = r. 

Note 1: The FLOOR(SQRT(x-10)) function can easily be modified to show that a function is unde- 
fined on other subsets of the real numbers. 

Note 2: All of the above functions PLOT satisfactorily in my DFW Version 4. 
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Francisco Cabo Garcia, Spain pcabo@esgueva . eco . uva . es 

We have a problem with DfW 4.06. 

If we ask the computer to calculate SIN(ASIN(x)), the answer is x but when we ask to calculate 
ASIN(SIN(x)), the outcome is: 


SIN (ASIN (x) ) = x 
ASIN(SIN(x)) = SIGN (COS (x) ) 


We would like to know whether it is possible to ask DERIVE to return x instead of the expression 
above? 


/ 

x 

1 1 

\ 

x - n« FLOOR 

1 


1 


V 

l n 

2 J 

/ 


Dr. Nigel Backhouse, Liverpool sx52@liverpool .ac.uk 

Francisco, 

The fact is that ASIN(SIN(x)) does not coincide with the function x. It is a sawtooth type of function, 
with discontinuities. Try to plot the output function. 


(I tried, and you can see how it looks like in DERIVE 6.10:) 


ASIN(SIN(x)) = SIGN (COS (x))- 


^jt-FLOOR^ 


1 

2 




Johan Vegter, Enschede, Netherlands cad0hartman . nl 

Hello DERIVERS, 

Em a relative new user of DERIVE and want your help with a (small) problem: 

I have two equations: 


(i) 


. 360 (L + jt) r R 

sin = L — 

AnR 2 


( 2 ) 


f = R 



f-1 

UJ 


2 


I am interested in the relation between / and x. R is unknown and L is a constant value. So, the idea is 
to solve this set of equations for R , and make one (large) equation with/ x and L (R is eliminated). 
Thanks in advance. 
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Terence Etchells, Liverpool T , A , Etchells@livjm . ac . uk 

Hi Johan, 

it is not possible to express f in terms of L and x (mathematics is the problem not Derive). However 
we can express x in terms of f and L, so you can find particular f s given L. I attach a MTH file that 
goes through steps. Hope this helps. 

Cheers, Terence 


#1: CaseMode := Sensitive 

Equation 1 and Equation 2 

#2: SIN 


r 360 -CL + x) L-F 

L 4-tt-R J 2 

-{■’-(tT) 


This is Terence’s procedure per- 
formed with DERIVE 6.10. 


#3: f = R 

Square both sides of #3 and then solve for R 


# 4 : 


H-f - (tAF ■ [* - ° 


2 2 "i 

2 4-R - L 


f 2 2 1 

2 4-R - L 

#5: SOLVE (R - f) = r R = 


r 2 2 

L + 4-f 

R = 

V. S-f 


Substitute for R in Equation 1 and then simplify 

■ R 

2 


#6: SIN 


r 360 -CL + x) h L-F 

L 4-tt-R J 2 


#7: SIN 


{ 720- f-x 720 -L -f 

+ 


2 2 2 2 
L tt-(L -t- 4-f ) ir-(L ir 4-f ) j 


2 2 
L - CL + 4-f ) 

16 -f 


Solve for x 




SOLUTIONS 


SIN 


( 720-f-x 720 -L -f 

+ 


2 2 2 2 
L jt-(L h r 4-f 3 rr-(L + 4-f 3 j 


2 2 
L - CL + 4-f 3 

16 -f 


# 9 : 


2 2 

jt-CL t 4-f 3 -ASIN 


2 2 h 

L-(L + 4-f 3 


16- 1 f I 


2 2 

TT- (L Hr 4-f J-ASIN 


720- | f| 

2 2 ^ 
L - CL + 4-f 3 


16- | f | 


720- f 


2 2 

2 2 2 2 rr-CL + 4-f J-ASIN 

tt -L - 720 -L -f + 4-rr -f 


2 2 2 2 
TT -L t 720-L-f t 4-rr -f 


720 -f 

2 2 
L - CL + 4-f ) 


16- f 


720 -f 


720- I f| 
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Hugh Porteous, England 


H.L.Porteous@shu.ac.uk 


If you integrate SIGN(SIN(x))*SIN(x) from 0 to 2 tc then you correctly get the answer 4. However, if 
you declare n to be an integer and integrate SIGN(SIN(x))*SIN(tzx) from 0 to In then you get the an- 
swer 0. Does this mean that DERIVE does not consider 1 to be an integer? 


Ray Girvan, England ray . girvan@zetnet .co.uk 

http : /_/www .users . zetnet . co . uk/rgirvan/ 

I'm probably missing something, but why would this imply DERIVE not considering 1 an integer? 
The answer does surprise me: with n otherwise undefined, I can't see why DERIVE returns 0, rather 
than a function of n. 


2-7 T 

/ SIGN (SIN (x) ) ■ SIN(x) dx = 4 

0 

n :e Integer 
2 -7T 

/ SIGN (SIN (x) ) ■ SIN(n ■ x) dx = 0 

0 

If you try the integral with n real, the answer comes out as fiji) .... 

n : c Real 
2n 

COS ( 2 • n • n ) 

SIGN (SIN (x) )• SIN (n-x) dx = 

0 n 

COS ( 2 • n • n ) 2-C0S(n-n) 1 

- + 

n n n 


2 • COS ( n • n ) 1 

+ — 

n n 


which gives the required answer for integer values of n. Particularly, f{n=l) = 4. 

This is the DfW 4 result. See below the DERIVE 6 result, 

n :e Real 


2 -7T 

J SIGN (SIN (x)) - SIN(n ■ x) dx = 
0 


COS(tt - n) ■ (C0S(2 -7T- n) - 1) 
n ■ (COS(tt- n) + 1) 


SUBST 


( COS (jt - n ) - (COS (2 « 7T ■ n ) - 1) 
n - (COS(tt - n) + 1) 


■p n ? 1 


COS(jt - n) - (C0S(2 -7T - n) - 1) 

lim = 4 

n+1 n-(COS(rr-n) + 1) 


The results are not identical! 
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n 

C0S(2 - 7T - n) - 2 ■ C0S(tt - n) + 1 CGS(tt. n) - (C0S(2 -tt- n) - 1) 2 2*(-l) 


n 


n ■ (COSCtt - n) + 1) n n 


Francisco Marcelo Fernandez marcelo0diquite . edu . ar 

I obtained the following results for the integral SIN(x)*SIN(tz x) from 0 to 2n: 

1) 0 when n is defined to be an integer 

2) 7i when n= 1 

3) SIN(2tt n)/(n A 2-1) when n is real. This expression yields the result n when n — > 1 


Johann Wiesenbauer, Vienna, Austria J. Wiesenbauer0tuwien . ac . at 

Hi all, 

Frankly, those results aren't that unexpected in my eyes. In the first place, assuming that n is any real 
number the exact result in DERIVE-notation would be 

IF(t 7 = 1, 7i, \f{n = -1, -71, SIN(2 n n)/((n + 1 )*(« - 1)))) 

as easily can be checked. As on many other occasions (try e.g. the simple example INT(x a t/,x)!), 
DERIVE gives only the general term and leaves the treatment of the special cases to you. This is a 
compromise I for my part can live with, because I can imagine that otherwise things might become 
fairly complicated if this result with its nested IF-statements is not a final one but a subterm in a much 
bigger term! 

If n is declared as integer then obviously DERIVE simply replaces the numerator SIN(2 7i n) of the 
general term by 0, wich leads to the value 0 for the fraction in general, but as before the cases n = 1 
and n = -1 are exceptions which aren't considered by DERIVE. (Interestingly enough, even Marcelo 
seems to have overlooked the latter case!) 

Thus in my opinion, unlike the original example by Hugh Porteous which is still a puzzle to me, the 
result given by DERIVE in this case is exactly what you can expect of a CAS. 

Cheers, Johann 


Lorenz Jaeneke, Switzerland ljaeneke0datacomm. ch 

The result of the considered integral INT(SIGN( ,0,27i) is indeed the expression given below, 

which leads for some integers to 


VECTOR 


f COS ( 2 • n • n ) 
v n 


2 • COS (n • n ) 1 

+ , n, 0, 10 

n n 


4 4 4 4 

?, 4, 0, , 0, , 0, , 0, , 0 

3 5 7 9 


The quotation mark tells us about the both vanishing numerator and denominator. So, we have to cal- 
culate the limit: 


r COS ( 2 • n • n ) 2 • COS (n-n) 1 N 

- + 


lim 

na>0 


n 


n 


0 
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All these results are clearly visible in the plot of the function f(n) = COS(2 tt7z)/?z - 2 COS(nn)ln + 1 In. 
Just as clearly as when we compare the plots of SIGN(SIN(x)) SIN(^zx) for even integers with those for 
odd integers (pay attention to the symmetry). So it's not an error it's a feature. 



iTiT^nY fT? Y F]:t Y FHt y f? Y ' 

t j— Rlgebra Calc Other FrgnlQ Clean Up 


■ j Q i. :=. i gm, £ i m, y j j ■ £ l nm ■ yjj a x. | n = i 


1 Jg 71 (sign(sinCx)) ■ sin(n ■ x))dx | n = 2 

4.S 

1 Jg 71 (sign(sinCx)) ■ sin(n ■ x))dx | n = 3 


4.00 


4.82e-14 


1.33 


■■■in<x>>ttsin<nttx> J .x J .0 J .2^n> I n=4 

hi H I H F;hD EKHCT FUNC 


Let’s compare earlier DERIVE results with recent ones: 


( C0S(2 -JT - n) - 2 ■ C0 S(tt - n) +1 ^ 

VECTOR , n, 0 P 10 


4 4 4 4 

? ? 4p , 0, f 0, , 0 P f 0 

3 5 7 9 


( C0S(jt - n) ■ (C0S(2 - tt - n) - 1) ^ 

, n, 0, 10 

n ■ (C0 S(tt- n) + 1) 

[ 7 , 7 , 0, 7 , 0, 7 , 0, 7 , 0 p ?, 0] 

COSi>.rO-(CQSC2.7T.rO - 1) 


VECTORl lim 
*k 


{ v 


n ■ (COSCtt- n) + 1 ) 


k , 0, 10 


4 4 4 4 

0 P 4 P 0 S , 0 , , 0 , , Op , 0 

3 5 7 9 


Harald Lang, Stockholm, Sweden lang0math . kth . se 

I have DERIVE 3.14 classic. When solving i SIN(w) - 3 COS(w) = 2 for w DERIVE produces the 
erroneous answer: 

JT f n ^ 

w — + LA LN(V 2 ) 

2 L 4 J 


Do later versions of DERIVE solve this equation correctly? 
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DNL± 

See DERIVE version 4 's answer: 

i-SIN(w) - 3 • COS ( w) = 2 
SOLVE (£• SIN (w) - 3 • COS ( w) = 2, w) 
3 • n £ • LN ( 2 ) 


3 • n 


4 2 4 

DERIVE 6 gives the same results. 


i . LN ( 2 ) 5 • n 

, w = 


i • LN ( 2 ) 


Steve Smith, Powell, Wyoming 

I am trying to plot Euler's Spiral: B(t) = 


r 

( 2 




r 

( 2 


COS 

n x 


dx 

+ i 

sin 

nx 


2 

J 

2 

J 






V 


J 


V 


smiths0mail . nwc . whecn . edu 


dx . 


When I try to plot the following: 

[LNT (COS (tt/2 x A 2 ) ,x,0 ,t) , INT(SLN(tt/2 x A 2),x,0,t)] 

DERIVE says ’’sorry, the highlighted expression cannot be plotted”. I assume it is bailing out because 
of the integrals, but I don't see why that should be a problem. I can compute a few of the individual x 
and y coordinates using vector expressions such as 

VECTOR (LNT (COS (tt/ 2 x A 2) ,x,0,t) ,t,0,10) 

VECTOR(LNT (SIN (tt/ 2 x A 2) ,x,0,t) ,t,0,10) 

But combining the results by hand for DERIVE to plot is time consuming. Is there a nice way for 
DERIVE to plot this type of function? 

Thanks, Steve 


Peer van de Sanden, Netherlands 

Hello, Steve, 

I have copied your vector to DfW 4.1 and tried to plot 
it. No problem at all. Then I tried DERIVE XM 3.05 
for DOS and tried to plot the vector. It was slow, but 
this version has also no problems with the plot. 


No problems to plot the parameter representation 
of the spiral with DERIVE 6. 


peer . mari j ke0gironet . nl 
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The Erlang C formula story 


Ray Girvan, England ray . girvan6zetnet .co.uk 

I've just offered (rashly, it turns out) to confirm, with Derive in exact mode, some results for the Er- 
lang C formula. 

E (m, u, rho) : =u A m/m ! / (u A m/m! + (1-rho) *SUM (u A k/k! , k, 0,m-l) ) 

My timings, using DFW on a 133MHz Pentium, are: 


E(320,3 19,0.9999) ... 8 seconds 
E(640, 639, 0.9999) ... 68 seconds 
E(1280, 1279,0.9999) ... 1098 seconds (ie 18 minutes) 
E(2500, 2499, 0.9999) ... 5359 seconds (ie 89 minutes) 


(0.406 with DERIVE 6.10) 
(2.75) 

(5.55) 

(15.9) 


I'm wondering how much can be trimmed off the time. Is there anyone with a faster machine who'd be 
prepared to try out some/all of the above for timing comparisons? Thanks in advance, Ray 


Thomas Lingefiard, Sweden 

My timings on a 266 MHz Pentium II Notebook (80 MB RAM; Win98, DfW 4.1 1) 

E(2500, 2499, 0.9999) ... 2464.4 seconds (i.e. 41 minutes 4 seconds) 

On a 450 MHz Pentium II Workstation (264 MB RAM; WinNT, DfW 4. 1 1) 
E(2500, 2499, 0.9999) ... 1375.6 seconds (almost 23 minutes) 


Dr N.B. Backhouse. Eneland 

Timings on a PC fitted with a 300 MHz AMD K/6-2 processor (64 MB RAM; Win98; DfW 4.1 1) 

E(320,... 1.9s 

E(640,... 15.5s 

E(1280,... 150.5s 

E(2500,... 1109.2s 

Impressive! 

Johann Wiesenbauer J . Wiesenbauer@tuwien . ac . at 

Hi all, 

Rather than using faster and faster Pentium PCs it may be worthwhile to look into Ray's formula itself 
which is certainly a far cry from being optimized. For example, the simple approach of rewriting it in 
the form 

E (m, u, p) : =u A m/ (u A m+m ! ( 1 — p ) LIM (TAYLOR (EXP (x) , x, 0 , m-1 ) , x, u) ) 

will already about half the computation time, but there should be even more time saving implementa- 
tions. How about a competition on that score? 


Cheers, Johann 
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Stefan Welke, Germany Spwelke6aol . com 

Hallo derivers, I can offer an improvement of Ray's Erlang C formula. This is my version: 

AUX(x , n) := n • ITERATE([1 , 1] + v_il.[l, v_i2/x] , v_, [1, 1] , n - l)/x 
E2(n, x, r) := 1/(1 + (1 - r)-(AUX(x, n))i2) 

I found the following execution times on my Pentiumll, 300Mhz: 

E(640, 639, 0.9999) : 24.3 seconds (Ray's version) 

E2(640, 639, 0.9999) : 6.9 seconds (My version) 

E(2500, 2499, 0.9999) (0.39 sec with DERIVE 6.10) 

This is about twice as fast as Johann's version utilizing the Taylor polynomials (5.34 sec with 
DERIVE 6.10). 

Cheers, Stefan 


Johann Wiesenbauer 


J . Wiesenbauer0 tuwien . ac . at 


Very nice, indeed! Although I would prefer a more concise representation such as 


E3(m, 


LI 



f 

f 

v_ 





2 


li + (1 - p)-m- 

ITERATE 

[1, 1] + V_ ■ 

i, — 

, v_, [1, 1], m - 1 



1 

LI 

Jj 


from an algorithmic point of view this could well be the very best that can be achieved. Isn't it con- 
solatory that in the first place it's the gray matter of the programmer that counts and not the number of 
MHz of his machine? 


Cheers, Johann 


From Josef Bohm to Ray Girvan: 

I have a very stupid question: For which purpose do you need Erlang C? I looked into several 

formula collections, but nowhere Erlang C is even mentioned? 

Ray Girvan, England ray . girvan6zetnet .co.uk 

Josef, 

This has turned into a very interesting problem (and thank you, everyone, for the improvements). 

I had never heard of it either. It's something to do with queuing theory in telecommunications, though 
I’m not sure what, as I don’t use it myself. The reason for needing to test this was to check a Mathcad 
worksheet I wrote. Someone in the Mathcad discussion group was having problems with this func- 
tion, since Mathcad's symbolic solver was falling over, and its numerical solver overflows at 10 A 307. 
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I tried recasting the Erlang C function as an approximation using the logarithm of the Stirling formula, 
so that all the power and factorial work could be done in the log domain to keep the magnitude safe. 

LOG_GAMMA ( n ) : =LOG (SQR (2*pi) ) -n+ (n-0.5) *LOG(n) +LOG(l+l/ (12*n) +1/ (288*n A 2) - 
139/ (51840*n A 3) -571/ (2488320*n A 4) +163879/ (209018880*n A 5) + 

52 4 6819/ (75246796800*n A 6) -534703531/ (902 961561 600*n A 7) - 
4483131259/ (86684309913600*n A 8) ) 

F (k,m, u) : =EXP (LOG_GAMMA (m+1 ) + (k-m) *LOG (u) -LOG_GAMMA ( k+1 ) ) 

ERLC (m, u, rho) : =1/ (1+ (1-rho) *SUM(F (k,m, u) , k, 0,m-l) ) 

It's messy and inefficient in Derive, but runs extremely fast in Mathcad, to an acceptable 13 or 14 dig- 
its of accuracy. For instance, 

ERLC(5000, 4999, 0.9999) = 0.991 149056003 13 ... 14 seconds! 

(I think this speed is due to the way Mathcad throws away small terms below the set tolerance level). 
But my problem then was to prove its accuracy against exact solutions. 

I have to admit, when I first saw the name Erlang C, it sounded like a dialect of C developed at Erlan- 
gen! I did a small web search; Erlang turned out to be a person, who developed loading formulae for 
telephone exchanges early this century. See: 


http : // http : //plus .maths . org/content/os/issue2/erlang/ index 

Ray 

End of the Erlang C formula story 
Comment: the URL given above is valid, Josef 


Subject Pythagorean Triples (Requested by T.P. Woods) bjaiien@ibm.net 


I note that all Pythagorean Triples can be commutatively expressed, as the sum of consecu- 
tive, odd, integers, two within the confines of the third e.g. 5 A 2 = 3 A 2 + 4 A 2 = 4 A 2 + 3 A 3 = 
= [1+3+5+7+9] = [1+3+5] + [7+9] = [1+3+5+7]+[9], by virtue of the commutative property of 
addition, and as correspondingly bracketed. 

Similarly 13 A 2 = 5 A 2 + 12 A 2 = 12 A 2 + 5 A 2 = [1+3+5+7+9+11+13+15+17+19+21+23+25] = 
= [1+3+5+7+9] + [11+13+15+17+19+21+23+25] = [1+3+5+7+9+11+13+15+17+19+21+23] + 
+ [25] as correspondingly bracketed, and so on for all Pythagorean Triples, and there is an 
infinite number of them. 

Surely, there must be a formal algebraic proof of this simple commutative property that is 
unique to squares, and never to cubes, or any other higher power? 

I have searched everywhere without success. Perhaps, your mathematics department may 
contribute. I'd appreciate your help. Thank you. T.P. Woods 
(From the cti-maths Newsgroup) 
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Do you know “Vampire Numbers”? 

Project: The Great Canadian World Wide Web Vampire Number Search 

Project Overview: My grade 12 computer science class and I are looking for about 1,000 participants 
who will run a small computer program overnight and email back a data file of results. We would like 
participants from all over the world. 

Details: A vampire number is a whimsical mathematical idea, introduced by computer graphics guru 
Clifford Pickover, in the June 1995 issue of Discover magazine. It makes note of the fact that when you 
multiply two numbers, all of the digits in those two numbers occasionally show up in the result. For ex- 
ample 21 x 87 = 1827, or 146 x 938 = 136948. The number 1827 is a four digit vampire number and 
136948 is a six digit vampire number. Looking for these numbers makes for an interesting computer 
programming project. I have several web pages that describe vampire numbers in detail. See 
http://grenvillecc.ca/facultv/ichilds and follow link(s) to vampire numbers. (Or you can go directly to 
http://grenvillecc.ca/faculty/ichilds/vamlil2.htm ) 

Comment: the URLs given above are not valid. If you are interested in Vampire Numbers then 
try the links given below, Josef 

http://en.wikipedia.org/wiki/V ampire number 
http://rosettacode.Org/wiki/V ampire number 

http://www.dreamincode.net/forums/topic/204318-c-program-to-generate-all-vampire-numbers-upto-a- 

given-number/ 


Riideger Baumann, Celle, Germany baumann-celle6t-online . de 

Dear Josef, 

The "Look'n Say"-Sequence {1,11,21,1211,111221,312211,....} seems not to be easily to program in 
DERIVE because you have to consider the order of the numbers. 

The sequence which I presented in DNL#33 {1,3,7,12,18,26,35,...} forms together with the sequence 
of its differences [2, 4, 5, 6, 8, 9,....} the set of the natural numbers. It is not difficult to program this se- 
quence in PASCAL. 

Another thing: Here is a DERIVE session belonging to Alfonso's ACDC7 and Leandro's and Javier's 
contribution in DNL#30: Diophantine Equations: 

#1: InputMode :z Word 

#2: El algoritmo de Euclides 

#3: EUCLIDES(a, b) :z ITERATE(IF(y = 0, [x, y, sO, si, tO, tl], [y, M0D(x, y), 

si, sO - FL00R(x, y)*sl, tl, tO - FL00R(x, y)-tl]), [x, y, sO, si, tO, 



tl], [a, 

, b, 1, 0, 0, 1]) 

#4: 

LET(u, v, 

w) := ITERATECu, v, w, 1) 


EUCLID(a, 

b) LET(T v , v , v 1, v, EUCLIDES(a, b)) 

#5: 


L 1 3 5j 
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#6: 

Eje 

implo de 

Tortosa 

& Santacruz 

#7: 

[a 

:= 29, b 

i r 

■--j 

n 

:= 7 SO] 

#8: 

EUCLID(a, 

b) = [1, 

-7, 12] 

#9: 

[d 

:= 1 s : 

1 1 
1 

-vl 

r+ 

1 1 

H- 

r-j 

i i 


#10: 
# 11 : 
#12: 
#13 : 


Solucion particular: 



Solucion general : 



j = p0 [-5460 , 


- u [17, -29] 


9360] 


#14: 

p0 + k ■ LI = [17. k - 5460, 9360 - 29-k] 


#15: 

2 

(F(k) CpO + k ■ u) ) = F(k) : 

2 

= 1130. k - 728520. k + 117421200 

#16: 

SOLVE (2260 -k - 728529, k) = 

r 728529 ^ 

k - 


ft. — 

l 2260 j 


#17: 

SOLVE (2260 -k - 728529, k) = 

(k = 322.3579646) 

#18: 

kO 322 



#19: 

ScHiicion special : 



#20: 

pO + kO-Li = [14, 22] 



#21: 

|29*14 + 22 ■ 17 = 780H-. :y. • 




Jonathan Devor, USA j devor6usa . net 

Hello all DERIVERS, 

I've been trying to solve the following set of non-linear equations: 

b 2 + b 'C + c 2 = x, a 2 + a -c + c 2 = y, a 2 + a -b + b 2 = z 

solving for a, b, c. 

I use DERIVE 3.12 for DOS, and as is, it gave up on these equations very quickly. Helping it along 
manually, I got a ’’solution not verified” or an unending calculation pause (I gave up after 20 hours on 
my 266MHz Pentium). Does any one out there have a substitution or other trick to get a solution for 
this, seemingly simple problem. 

Thanx. 

Johann Wiesenbauer J . Wiesenbauer6tuwien . ac . at 

In the enclosed file I have tried to sketch an ad hoc solution (without Groebner bases and all that 
stuff!) If you feel like trying it out yourself, you are strongly advised to print it out and input or copy 
all the necessary lines into a new (!) Derive-file. (Note that double assignments must always be simpli- 
fied subsequently, otherwise they won't work properly!) 
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I have left out some special cases and details, but would be glad to provide you with further informa- 
tion if needed. One concluding question: Is this ’’seemingly simple” problem, as Jonathan put it, really 
simple or not? 

I leave it to you to answer this question. 

Cheers, Johann 

x = V6-V(V3-V(- a A 4 + 2 • a A 3 • (b + 2-c) - a A 2-(b A 2 + 2-b-c + 6-c A 2) + 2-a-c-(b A 2 - 
b-c + 2-c A 2) - c A 2-(b A 2 - 2-b-c + c A 2))-ABS(b - c) + 2-a A 3 + a A 2-(b - 5-c) + 
a-(b A 2 - 3-b-c + 4-c A 2) + 2-b A 3 - 5-b A 2-c + 4-b-c A 2 - c A 3)/(6-./(a A 2 - a-(b + 
c) + b A 2 - b-c + c a 2)) 

(This one solution for x- there are four. The solutions for y and z are obtained by substitu- 
tion. Johann changed the unknowns a,b,c to x,y,z and the parameters x,y,z to a,b,c. See 
NONLIN.MTH.) 

Ray Girvan, England ray . girvan6 zetnet .co.uk 

This is not as rigorous as Johann's solution, but a brute force approach seems to work (in DfW) 


#1: 

b A 

2 

+ 

b*c 

+ 

c A 

2 

- X 

#2: 

a A 

2 

+ 

a*c 

+ 

c A 

2 

- Y 

#3: 

a A 

2 

+ 

a*b 

+ 

b A 

2 

- z 


Solve #1 for b gives two solutions of form b = f(c,x) and solve #3 for b gives two solutions of form 
b = f(a,z). Pair all combinations of these solutions to give four equations of the form f(a,c,x,z) = 0. 
Solve each of these for c, giving eight solutions c = f(a,x,z). 

At this point, plug in trial values .... which shows four are duplicates, which you can reject. Then plug 
these c results into equation #2, simplify, and solve for a. This gives sixteen solutions, which again 
contain duplicates. They look like, for instance: 

a = SQR ( 6) *SQR (-SQR (3) *SQR (-x A 2+2*x* (y+z) -y A 2+z* (2*y-z) ) 

*ABS (x A 2-x* (y+z)+y*z)-x A 3+4*x A 2* (y+z) -x* (5*y A 2+3*y*z+5*z A 2) 
+2*y A 3+y A 2*z+y*z A 2+2*z A 3) / ( 6*SQR (x A 2-x* (y+z ) +y A 2-y*z+z A 2 ) ) 

And these can then be substituted back into 1#, #2, #3 to get expressions for b and c. 

Ray 

Johann Wiesenbauer J. Wiesenbauer6tuwien , ac . at 

Hi all, 

[Unfortunately there had been a lot of mistakes in my last email. For your convenience, what follows 
is a (hopefully!) correct version. I'm very sorry!] 

Since I have some time left right now, some overdue comments on Jonathan's "simple” set of nonlin- 
ear equations might be in order. In the first place it shouldn't be too difficult to come to the agreement 
in the community of mathematicians that it is more appropriate to use the following notations: 

x A 2 + xy + y A 2 = a 

x A 2 + xz + z A 2 = b 

y A 2 + yz + z A 2 = c 

(Otherwise even DERIVE gets confused in this topsy-turvy world, where a,b,c are variables and x,y,z 
constants!) Anybody out there who objects to it? 

Here we go! 
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This is actually a quadratic equation for w A 2 (I hear Jonathan moan in the background!), so you can 
solve it even without DERIVE. One solution is 

wl := sqrt ( (sqrt (3 ( (a+b+c) A 2-2 (a A 2+b A 2+c A 2 ) ) ) +a+b+c) / 2) 

and the other three you get by changing the signs of the two square roots independently. The corre- 
sponding values of u and v are now given by 

u = (a-c)/w, v = (a-b)/w 
and x,y,z are finally obtained by the formulas above. 


One more word to Jonathan's linear equation 

(a-b)x+ (c-a)y + (b-c)z = 0 

which I have also found myself. I can understand that he enjoyed this beautiful equation, but it's more 
interesting from a philosophical point of view (why are all solutions - typically 4 as we have seen - on 
this plane and what does it mean from a geometric point of view?) and doesn't provide that extraordi- 
nary shortcut in my eyes. (Furthermore, it doesn't exist in the case a = b = c, as Ray has noticed!) Look 
at my previous DERIVE-solution again and you'll see that I arrived very fast (and without any manual 
computations!) at a fourth equation of the form (#17) 

x A 2 (a - 2b + c) - xy(2a - b - c) + y A 2 (a + b - 2c) = - (b - c) A 2 

This is what you can also expect when using this linear equation. But the actual struggle began there- 
after! (Well, with the benefit of hindsight maybe I should have covered this up by saying "the rest is 
plain sailing"!) 


Cheers, Johann 

DERIVE 6 has no problems solving this nonlinear system using its built-in Groebner bases. 
The two lines of code given below produce all solutions within 17 seconds, Josef 

2 2 2 2 2 2 
Leql := b + b ■ c + c = x, eq2 := a + a-c + c = y, eq3 := a + a«b + b = zj 

SOLUTIONS (eql a eq2 a eq3 , [a, b, c] ) 


Arnaldo Struzberg, Rio de Janeiro, Brasil aesse6marlin . com.br 

How can I obtain with DERIVE the cartesian graph of a binary relation? 

Examples: S = {(x,y) e RxR| x A 2 y < y A 2} 

S = {(x,y) e R x R | x A 2 + y A 2 < 4} 


Thanks 

DNL: no j o . boehm6pgv . at 

What about this? Here is an attempt to graph your relation as a collection of points. Approxi- 
mate and plot. Josef 

There are rumors, that DfW5 will be able to plot the graph of a relation directly! In the mean- 
while I could test DfW5 and I can confirm, it does. 
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REL(rm, xa , xe , dx , ya , ye, dy) := VECTOR (VECTOR ( IF ( rm , [x, y] , ?, ? 

), 

x, xa, xe, dx) , y, ya , ye, dy) 


rm : 

: = X > 2 

A X 

1 = 

1 y 1 


REL i 

(rm, - 1 0 

, 10, 

1 , 

-10, 

10 , 1) 

rm2 

:= x > 

-2 A 

1 x 1 

< 1 y 1 


REL i 

( rm2 ? -5 

, 5, 

0.2 

, -5, 

5, 0.2 


REL(y > | x — l| a y < |4 — x|, -10, 10, 0.5, -10 
REL (MOD ( y , x) = 2, -10, 10, 1, -10, 10, l) 

2 2 

REL ( x + y <4, -3, 3, 0.1, -3, 3, 0.1) 

2 2 

REL(x • y < y , -5, 5, 0.1, -5, 5, 0.1) 


Above you can see the DERIVE 4 procedure but ... 


, 10, 0.5) 



... rumours were right. You can now plot the binary relations directly by simply entering the 
relation and the plotting it, Josef 



See how Tl-Nspire reacts on some User Forum requests from 1999 on the next 
pages. 
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arcsin(sin(x)) does not simplify further. 
Johan Vegter’s system can be solved! 



On the next page you will find that 

(1) Ray Girvan’s integral can not be given in general form, 

(2) Harald Lang’s equation can not be solved exactly, only numerically, 

(3) Peer van de Sanden’s spiral can be plotted (needs some time). 
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solvefr 

sin (tv) - 3 ■ cos (tv) = 2, tv) 

false 

c Solve ( 

i ■ sin (tv) - 3 ■ cos (tv) = 2, tv) 3 ■ cos (tv) + - r sin (t 

v)+2=0 

cSolvef 

i ■ sin (tv) - 3 ■ cos (tv) = 2, tv) i jt 



3 - cos (tv) +e ^ ■ sin(T 

+ 

to 

II 

0 

a c Solve ( i ■ sin (tf) - 3 ■ cos (tv) = 2 } w) 


w=2. 356 19449-0. 34657359- i 


3-jt i- ln(2) 2.35619449-0.34657359 - i 

4 2 

«' c Solve works only by approximation (ct!r+ Enter), no exact* 
B 



p 24 


TI-Nspire ’’ s Answers to the User Forum 


D-N-L#34 


The Erlang C formula can be defined and works in a very restricted amount. 



I didn’t try solving the nonlinear system from page 17. 

But I started developing a package for solving linear difference equations of first and second 
order. DERIVE provides a utility file RecurrenceEquations.mth for finding exact solutions of 
first and second order difference equations. 

It started very promising. Then I discovered an expression ATAN(V(-d),-p) which is explained 
in DERIVE’s help file as follows: “ATAN(y, x) is the angle of the point (x, y) in the x-y plane 
measured counterclockwise from the positive x-axis.” This was no problem to realize with 
Nspire. Next difficulty was the reason to give up - at least for me: 

h-Z(r-LIM([- ELEMENTS, 2), ELEMENT(h, 1)]/DET([h, LIM(h, x, x + 1)]), x, x + 1), x) 

This sum(u,x) gives the “antidifference” F(x) of u with respect to x. F(x+1 ) - F(x) = u 

3 2 

2 xx 7 ■ x 

I(x + 1) = + 

x 3 2 6 

-X - 1 -X 

-x - 2 -x 5-3 .(-1) 

110-3 -(-I) = 

x 2 


3 2 

(x +1) (x + 1) 7 ■ (x + 1) 

- + 


3 2 

- Cx + 1) - 1 


3 2 

x x 7 ■ x 

- + 


2 

x + 1 


5-3 


■(- 1 ) 


Cx + 1) 


3 2 6 J 

-x - 1 -x 

5-3 .(-1) -x - 2 -x 

= 10-3 -C-l) 


I was not able to find the antidifference (which is the indefinite sum) in my first attempt with 
TI-Nspire. Any advice would be highly appreciated. See more on page 57. 
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WHAT FILOU HAS BEHIND HIS HEAD 

Jean-Jacques Dahan, Toulouse, France 

Jean-Jacques DAHAN is a Mathematics Teacher in Toulouse (France) and the author of a 
book called "An introduction to geometry with the TI-92" in which he has created a small 
character baptized "Filou" which he uses to show how one can approach the module of ge- 
ometry using TI-92 in a recreational way. 

You are going to attend a conversation between Jean-Jacques and Filou where Jean- 
Jacques is going to try to show to his little friend that even the most harmless of configura- 
tions can hide subtle results. It will perhaps convince the reader of the power of the relation- 
ship between Computer Algebra Systems and Geometry using the same instrument. 


WHAT FILOU 

HAS BEHIND HIS HEAD 



STORIES ABOUT PARABOLAS 


FILOU: Tell me, Jean-Jacques, what's the meaning of this title? How could I have something 
behind my head without even knowing it myself? 

JJD You forget that I created you and I am no doubt in a better position than you to know 

that you are concealing a magnificent parabola. 

FILOU: Really! a parabola, this parabola, I don't feel it, , neither in front of my head, nor be- 
hind it. You must explain this to me, there is something in what you say but I don’t 
know what and I don’t understand. 

JJD Of course! I’ll explain it to you as well to our readers who I feel, are spying in our 

conversation with certain curiosity. I am going to tell you a mathematics story in 
three parts and you will see that your head has a very close relationship with a cer- 
tain parabola ! 
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FIRST EPISODE: 

ALL POLYNOMIALS OF THE SECOND DEGREE LEAD TO A PARABOLA 

1/ DEFINITION: let's call a parabola any curve of the euclidian plane that is the representa- 
tive curve of a function of the type x -> ax 2 in the cartesian coordinate system of this plane 
with a^O 

21 CLASSIC PROPERTIES: Theorem 1 : the representative curve of a function of the type 
x^ax 2 +bx+c with a^O in a cartesian coordinate system (O, i, j) of the plane is a parabola. 

Proof: 

It suffices to seek the equation of this curve in the cartesian coordinate system (I, i, j) , where 

i_ 

I is its extremal point of the abscissa - — .. One finds Y=aX 2 . 

Z.3 

Theorem 2: a parabola is the set of the points M of the plane such that MF=MH where F is a 
point on the plane called the focus of the parabola and H is the orthogonal projection of M on 
a line D called the directrix of this parabola. If y=ax 2 is the equation of the parabola in a car- 

1 

tesian coordinate system, then F has for its coordinates ( — , 0 ) and D has for its equation 



Proof: it is an elementary norm calculation that can be made very easily with the TI-92. 

3/ A MORE ORIGINAL PROPERTY: 

Constructions, manipulations and conjecture 

Let's draw in a cartesian coordinate system 
(O.i^jJ, the curve of the equation y^ai.^) 2 
with, for instance a-i=1/2. It seems that this 
curve has an axis of orthogonal symmetry 
parallel to the y axis of our coordinate sys- 
tem. Moreover, one can conjecture that this 
curve is a parabola as has been defined 
above. 


Note that we detect this symmetry by tracing 
a parallel line to Oyi until it matches up. We 
then draw the symmetrical point mi' from the 
generic point ml of our curve with respect to 
this line. When we move mi on our curve, 
mi' also seems to remain on the same 
curve. 
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Proof: 

Let (P) be the curve of equation y^ai.^) 2 in 
the cartesian coordinate system (O.i^jJ; let's 
find its equation in the direct coordinate sys- 
tem (0,i 3 ,j 3 ). For that, let's use the direct 
coordinate system (0,i 2 ,j 2 ). 

Note that (i 1 ,j 1 ) = a. . 

From the vectorial equalities: L, =i 2 and j., =cosa i 2 +sina j 2 and i 3 =sina i 2 -cosa j 2 and 
j 3 =cosa i 2 +sina- j 2 one can easily deduce the formulas : 



x 2 =Xi+y r cos a 
y 2 = yi-sin a 


Then : x 2 =sin ax 3 +cos a y 3 

y 2 =-cos a-x 3 +sina-y 3 


Then : 

1 j cosa . 

Xi — T — ,x 3 and yi : — .x 3 +y 3 

sina sina 

All calculations made, (P) has an equation in (0,i 3 ,j 3 ) : 

y 3 = 3l .(x 3 ) 2 + C0Sa (x 3 ) , so that (P) is indeed a parabola. 
sin 2 a sina 

NB: if I is the maximum point of this parabola, I ( x 3 (l) = si3 ( 2a ) , y 3 (l) = ), the former 

4ai 4ai 

has an equation of y 4 = a \ .(x 4 ) 2 in the coordinate system (l,i 4 ,j 4 ) where i 4 =i 3 and j 4 = j 3 

sin 2 a 

will have as its focus: 

F( x 4 (F)=0 , y 4 (F)= ^ ) or F( Xl (F) = , Yl (F) = -±~) 

And as directrix, the line (D) has the equations: 

V4=^ or y,=(. cosa ). x , 1 l- . 

Case: taking any initial coordinate system (0,1,]) with (!;]) = a 

Let us use in this coordinate system, a curve of the equation y=a.x 2 with a*0. 

An equation of this curve in the coordinate system. (O.i^jJ where: i 1 = jp r jj = anc * Ji = jpjj ' s 

yi =a1.(x1) 2 with al = 

Hill 2 


given by the formulas: x = jpjj anc ' Y = jpjj- It is therefore a parabola and we have proved: 
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Theorem 3: the curve of a function of the type y=a.x 2 in any coordinate system (0,i,j) of the 
plane with a^O and (i; j) = a is a parabola whose symmetry axis is a parallel line to the y axis 
of the coordinate system , this axis has for its equation : x = - c 5 >Sa , the focus F has for its 

23i 

rosr/ 1 

coordinates (x(F) =- ,~ , y(F) = - — ) and the director (D) has for its equation y=(-cosa).x- 

28i 48i 

1 3* Hi!! 

- — , where we have previously proposed al = i J . 

4a-! ||i|| 2 

End of the first episode 

FILOU: I think I’ve been quite patient up to this point, but where is the parabola that I’m sup- 
posed to be hiding ? 

JJD I’m getting to it; you should be more patient. This episode was just a necessary in- 

troduction. 

FILOU: Right! Well I guess I’ll have to summon up some patience. I’m all ears! 

SECOND EPISODE: 

A PARABOLA, FOUR TANGENTS AND PROPORTIONALITY 

Constructions, manipulations and conjecture 

Let the parabola of the equation y=-^.(3.^/3x 2 +6.x+^/3) be which is defined in the functions' 

editor under the name yl . This curve is then made to appear on one of the geom-etry files 
which I have named "courbyl" thus allowing it to be worked in on a geometrical way: note 
that at the same time as the curve yl appears, we get a tangent at a generic point that can 
be moved by using the cursor on the point. Note that we draw the following points on our 
figure : Oi(0,0), 0 2 (-1,1), Ai(1,0), A 2 (0,1+-\/3) and that the four lines (0i0 2 ), (AiA 2 ), (OiAi) 
and (0 2 A 2 ) are tangents to this parabola. 

The line drawn in bold is any tangent to the 
parabola cutting (OiA^ at IVh and (0 2 A 2 ) at M 2 . 
Measurements are made to the machine at the 
distances OiMi, OiA^ 0 2 M 2 and 0 2 A 2 . 


The ratios t^x 1 an d 9 2 ^ 2 are then displayed. It can be seen on the display that these ratios 

OiAi O2A2 

are equal and remain equal when one modifies the tangent used. The demonstration of this 
result, in this particular case, is relatively simple. 
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Equation of the tangent to our curve at Abscissa of the point 

the generic point M 0 (x 0 ,yo) 



Abscissa of the point Mi, which is con- Simple expression of this abscissa, and 

sidered as the intersection of the tan- expression of the abscissa of the vector 

gent at the point M 0 and the line 0 2 A 2 C) 2 ivf 2 which also proves that: 

whose equation is o» = 

y=V3x+1+V3 6 

The data of the four chosen tangents initially allows one to deduce, at least in this particular 
case that these are all lines (M-|M 2 ) such that Mi e (O 1 A 1 ), M 2 e(0 2 A 2 ) and = ^?iMi yet 

O2A2 OlAl 

respecting the relationship or non simultaneous relationship of Mi and M 2 to the half-lines 
[O1A1) and [ 0 2 A 2 ). 

One could re-construct the same experiment with four other tangents and end up with the 
same results. This results could also be obtained with the CAS module on the calculator. 

End of the second episode 

FILOU: Jean-Jacques, everything you’ve showed us is all very well, but I still don't see the 
relationship between my head and the parabola whose tangents you calculated with 
such skill. 

JJD: My dear Filou, If I may say so, you are suffering from a severe case of myopia. Look 

carefully at the figure again and tell me what you think of it. 









p 30 J.-J. Dahan: What Filou has behind his Head d-n-l#34 



FILOU: Goodness gracious! 

JJD And that’s not all! 

FILOU: With all due respect, may I point out that although this looks similar to my head, it is 
not my correct head. 

JJD Don't worry, be patient and wait until last episode. 


THIRD EPISODE: 

FILOU ! YOUR HEAD FOR A PARABOLA ! 


JJD. ..If youdon’t mind, I am going to use only the contours of your head ; your mouth and 
eye won't be useful. 


Constructions, manipulations and conjecture 

The point IV^ can move on the line (O1A1) 
when one pulls on the point positioned on the 
segment situated in the lower part of the 
screen (note that this segment is included in 
the axis of the abscissas of the hidden coor- 
dinate system). On the other hand M 2 is lo- 
cated on ( 0 2 A 2 ) so that: 0 2 M 2 = . 

0-|A.| 

The question that I put forward is the following: does there exist a curve which has the tan- 
gents MiM 2 when Mi moves along the line (O1A1). The calculator is going to answer in an 
experimental way using the option "Locus" obtained with in menu F 4 (one looks for the locus 
of Mi when the point situated on the lower segment covers this lower segment). We will 
therefore only obtain a small portion of the locus sought but this will suffice. 

It seems to be that one gets something that 
looks like a parabola. We must now prove it. 
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Proof: 



Choosing the cartesian coordinate system 
(0,i,j) as the one drawn on the figure on the 
left. The hypothesis of proportionality is 
translated into : 

0 1 M 1 =t-0 1 A 1 and 0 2 M 2 =t-0 2 A 2 


As M 1 M 2 • (a-t,p+t), an equation of the line (Mi M 2 ) can be written as: 
d(t) = 0 with d(t) =(x-t). (p+t)-y.( a-t ). 

One will obtain an equation of the envelope sought by eliminating t in the system: 
d(t) = 0 and d'(t) = 0, which we will do with the TI-92 : 



f ; ; -y 


f FE 1 

f s :: "I 


;! } i C 


PrgnlO 

L =!■■! :i ' -i: -_ I! .. 


■(x - t) (J6 + i) - y (oc - t) ■+ 


1 


1 solce( -2 ■ t + x + y - JS = G n i) t = 

, da) | t= 2i±pl 


■2-i + x + y- JS 
x + y - J6 


fil 


jZ r?-fr + R i-u 


d<t > I t=<x+y-B>/2 


fFiT^DT F£t 'I 
t f— Algebra 

f F3 t ■] 

Calc 

r fht 1 

Other 

f FE 1 

PrgnlO 

Ffi 

Clear a-z... 



#**(4*4)*#- 

■<* 4 )*#- 


(2-ffl + J8)-y JB 2 
2 4 


2 (2ct + JB)y | JB 2 


^-|x = z-y-JB 
z 2 

(-a-B)-y+4- 


/4-<2*«+ 0 >*y/2+ 0 A 2/4 1 x=h 


fihd e h ■: t 


FUNC lE/HO 


These calculations show, after the change of variable x=z-y-p, which corresponds to a 
change of coordinate system, that our envelope has for its equation in a suitable coordinate 
system: 

(-a-p).y+z 2 /4=0. It is therefore, in the general case, a parabola according to the first episode; 
the new coordinate system has for its origin O 3 (-p,0) and as its basis vector i and -i + j ; its 
symmetry axis is directed by -i + j . 

The particular case (-a-p)=0 corresponds to the case where (0i0 2 ) I I (AiA 2 ). 

1 

Finally, the equation which has allowed us to conclude, can be written: y = .z 2 

4.(a+P) 


Theorem 4: given the head of Filou... no. ..on a serious note, given a quadrilateral whose 
sides are not two by two parallel, the envelope of lines (M.,M 2 ) such as defined above, is a 
parabola tangent to the four sides of this quadrilateral. 

NB: if you want to know how I got this result using TI-92, you can read the brochure of the 
IREM of Toulouse that deals with this subject (among others things) entitled « De la propor- 
tionnalite aux coniques en passant par les lignes de niveaux avec la TI-92 » (IREM 1998). 


End of the third episode. 
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JJD: 

FILOU: 

JJD 


So, are you satisfied? 

I must confess that I am surprised to discover that my head like any other head 
could conceal such complicated calculaions, considering the fact that I’d alwayshad 
the impression that when looked at, it was after all my head. 

.We must now say to our readers that this article gives them the mathematical expla- 
nation of the pretty drawing which is on the cover of the book that I have published 
"INTRODUCTION A LA GEOMETRIE AVEC LA TI-92" (ELLIPSE Publishings) and 
which can be seen onthe front page of this article. 


Jean-Jacques DAHAN 
Mathematics Teacher. 
IREM of TOULOUSE * 


* Institute of Research for the Teaching of Mathematics. 


Many thanks to Henri LEMBERG who translated this article into English and above all to 
Caroline RENNIE who made the French to English translation comprehensible for Eng- 
lish speakers. 


ON THE SOLUTION OF A LINEAR DIFFERENTIAL 
EQUATION OF THE ORDER n WITH CONSTANT 

COEFFICIENTS (2) 

Emilio Defez Candel , Vicente Soler Basauri , Departamento de Matematica Aplicada 
Universidad Politecnica de Valencia, P. O. BOX 22.012. Valencia. Spain 

Example 3.3 Let us solve the following differential equation with DERIVE : 

/"(x)-5/'(x)+17j'(x)-13j(x) = 0 . (14) 

In this case vector v = [l, - 5, 17, - 13] , and we have: 

DIFERENCIAL_H0M0GENEA_C0EFICIENTES_C0NSTANTES([1, -5, 17, -13], x) 

2 ■ x 2 ■ x x 

e ■C0S(3-x)*d + e ■SIN(3'x)-r + e -c 

2,1 2,1 1,1 

2 ■ x 2 ■ x x 

Y(x) e ■C0S(3-x)-d + e *5IN(3-x)*r + e -c 

2,1 2,1 1,1 

Y"’(x) - 5-Y"(x) + 17-Y'Cx) - 13-YCx) = 0 


We can immediately check that the solution obtained is a solution of( 14). 
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Example 3.4 Let us solve the following differential equation with DERIVE - and check the solution: 

y (4) (x)-6/"(x) + 12/'(x)-8/(x) = 0 (15) 

In this case vector v = [l, - 6, 12, - 8, 0] , and we have: 


DIFERENCIAL_HI0MGGENEA_tQEFItIENTES_CQNSTANTES([l, -6, 12, -3, 0], x) 

2 2 ■ x 2 ■ x 2 ■ x 

x -e -c + x-e -c + e *c + c 

2,3 2,2 2,1 1,1 

2 2 ■ x 2 - x 2 ■ x 

Y(x) :=x-e -c + x-e -c + e *c +c 

2,3 2,2 2,1 1,1 

Y pppp fx) - 6-Y rpr (x) + 12-Y rr Cx) - S-Y'fx) = 0 


Example 3.5 Let us solve the following differential equation with DERIVE: 

y (4) (x) + 4 y"'(x) + 8 y"(x) + 8 /(jc) + 4y(x) = 0 (16) 

In this case vector v = [1,4, 8, 8, 4] , and therefore we have: 

DIFERENCIAL_HQMDGENEA_CQEFICIENTES_CQN5TANTES([1, 4, 3, 3, 4], x) 

-X -X -X -X 

x-e ■ COS(x) - d + e ■ COS(x) ■ d + x-e - SIN(x) ■ r + e ■ SIN(x) ■ r 

1,2 1,1 1,2 1,1 

-X -X -X -X 

Y(x) := x-e ■ COS(x) ■ d + e -COSCxJ-d + x-e -5INCx)-r + e -5IN(x)-r 

1,2 1,1 1,2 1,1 

Y prpr (x) + 4-Y rrp fx) + S-Y pp fx) + 3 ■ Y p (x) + 4-Yfx) = 0 
Example 3.6 Let us solve the following differential equation with DERIVE: 

y (4) (x) + 2 y"'(x) + 3 y "(x) + 2y'(x) + y (x) = 0 (17) 

In this case vector v = [1,2, 3, 2,1], and therefore we have: 

DIFERENCIAL _HOMOGENEA_CQEFICIENTES_CONSTANTES ( [ 1 , 2, 3, 2, 1] , x) 

- x/2 r ^3 ■ x ^ - x/2 C ^3 -x - x/2 f ^/3-x ^ - x/2 f ^/3-x 

x-e ■ d ■ COS + e ■ d ■ COS + x-e ■ r ■ SIN + e ■ r ■ SIN 

1,2 L 2 } 1,1 l 2 J 1,2 L 2 J 1,1 L 2 

Factor the above expression! 



Y pr,r Cx) + 2-Y ,rp (x) + 3*Y P P (x) + 2-Y P (x) + Y(x) = 0 
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3.1 Particular solutions. 


We can use the following function 

DIFERENCIAL_HOMOGENEA_COEFICIENTES_CONSTANTES(v,x ), 
in order to calculate particular solutions of differential equation (6) that verify 

y(x Q ) = e 1 , y'(x Q ) = e 2 , j y n ~ 1 \x Q ) = e n _ x 

For that purpose we will proceed like in the following example. 


Example 3.7 Let us consider the differential equation already solved in example 3.6 

y (4) (x) + 2y”'(x) + 3 y”(x) + 2y'(x) + y(x) = 0 . (18) 

In this case vector v = [1,2, 3, 2,1] and the solution written in a compact way is given above. 

Let us suppose that we are looking for the particular solution of (18) which verifies 

[y(o)=i,y(o)=2,y(o)=3,y"(o)=i]. (19) 

In this case, we calculate the successive derivatives of the solution function in the point considered, 
thus obtaining: 


[[YCO), Y r C0), Y"(Q), Y r r r (0)] ] ’ 


d 

1,1 


d = 1 

1,1 

d f3 ■ r 

1,1 1,1 
d - + 


d J'i-r 

1,1 1.1 

1,2 2 2 


d - + = 2 

1,2 2 2 

2-d + d - 2 * f3 * r + ^3 ■ r 

1,2 1,1 1,2 1,1 


2-d + d - 2 ■ J3 ■ r + J3 ■ r 

1,2 1,1 1,2 1,1 

2 


- = 3 

2 

3-d - 2-d + 3-73. r 

1.2 1,1 1,2 

2 


3-d - 2-d + 3 ■ 73 ■ r 

1.2 1,1 1,2 

= 1 

2 


We now proceed to equal each one of the results obtained to the conditions given in (19) in order to 
obtain an equation system (right matrix) which can be solved (using a substitution!): 


"dll = d, dl2 = e, rll = r, r!2 = s" 


SOLVE 



= 1, e 


d f3 • r 

+ = 2 , 

2 2 


2-e + d - 2*f3-s + f3-r 
2 


= 3 , - 


3-e - 2-d + 3-f3-s 


-6 a r - 


- 1 , [d, e, r, 
17- J3 


s] 


d = 1 a e = 


3 


a s = 2-f3 
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# 1 : 


# 2 : 


InputMode := Word 
r := rO- (1 - c- t), <p := 


c-t 


1 - c-t . 


#3 : v_ := 


#4: v_ := 


d d 

vr := — r , vtp := r ■ — cp 
dt dt 


c- rO 


vr := - c ■ rO , vcp := 


1 - c-t 



r d y 


2 

"d 1 

2 d 

#5 : a_ := 

ar := 1 — 1 r - r ■ 

— $ 

, atp := r* 

— 

tp + 2-vr- — ?p 


LdtJ 

L.dt j 


tdtj 

dt . 


#6: 3L := 


2 

c -rO 


ar := 


(c-t - 1) 


asp := 0 


(Comment: For further work it is necessary to simplify #3 and #5. This returns the velocity 
components in #4 and the acceleration components in #6.) 

To evaluate then the magnitudes of velocity and acceleration by entering and simplifying: 

#7: [ | [vr, v*p] | p | [ar s acp] | ] 


#8: 


Solution #8 will have simpler form if we declare the domains for parameters c and ro. We consider 
two cases: ro > 0 and c > 0 and ro > 0 and c < 0. 

For the first case: 

#9: [rO :e Real (0, *0, c :e Real (0, ™)] 




2 

2 2 

7(c *t - 2-c-t +2)- 

c-r0 

c ■ | r0| 

c-t - 1 

p 

3 


_ 


let - l| 


#10: 


2 2 

orQ-^tc -t - 2-c-t + 2) 

let - n 


2 

c -r0 


|ct - 1| 


Simplification of #8 returns #10. 
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Notes: 

We know that r component of the polar coordinates cannot be negative, hence from the content of the 
problem we have (1 - c t) > 0, so \c t - 1 1 = - (c t - 1). Using this result we can further simplify result 
#9 to the form: 

(Highlight | c • t - 1 | and substitute / simplify in one step - pressing CTRL+ENTER in DERIVE for DOS) 
(Comment: In DERIVE 6 you have to work via Simplify > Subexpression Substitution ....) 


# 11 : 


2 2 

c - r0-</(c ■ t — 2-c - 1 + 2) 
1 - c * t 


2 

c -r 0 


Cl - c-t) 

The same result is valid for c < 0 (do check it), but in this case the motion is unlimited in time. 


II.9 Trajectory of a moving point is given by the equations (in polar coordinates): 

r = Ae k °, O = Bt, (A, 5 , k - constants) 

Show that the angle between velocity and acceleration does not depend on time. 


Solution : The time independence of the angle implies the time independence of the trigonometric 
function of this angle (for example cos-function). So we will investigate the derivative with respect to 


time of the cos-function of the angle between velocity and acceleration: 

d 


r _ _ a 
v a 



_ 

_ 


V 

V 

a 

J 


d t 


First we enter formulae describing trajectory of movement, velocity and acceleration (some of them 
have already been entered in the previous problem - #2 and #3): 

[r := A-e A (k-t) , <p := B-t] 

Now entering and simplifying the derivative (#14 and #15 are results of simplifying #2 and #3): 


#12: 

CaseMode := Sensitive 

#13: 

[r 

:= i 

k.t 1 

f 0 := B ■ t J 

#14: 

v_ 

■A 

r k-t k-ti 

Lvr:=A.k.e , := A-B-e J 

#15: 

a_ 

A 

1" k-t 2 2 k«tl 

Lar := A-e ■ (k - B ), a<]3 := 2-A-B.k.e J 

#16: 

d 


[vr, v*p]"[ar, acp] 

- n 

dt 

i — i 
< 

“3 

1 1 

1 1 

1 

s 
1 1 

1 

c 


The result #16 proves the time independence of the investigated angle. 

Exercise: Try to solve the above problem using the cross product of two vectors. 
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11.13 Given are position vectors of three moving points: 

r x = cti + b y , r 2 = cti + bj + ek , r 3 = {ct + R cos(gH)) i + {b + R sin(cof)) j 
where c, b, e, R- constants, t - time. Prove that these points belong to the rigid body. 


Solution : The distance between any two points of the rigid body should not depend on time. Hence the 
equations: 


d 

dt 



= 0 


Let us enter the above vectors and simplify the appropriate derivatives. Result #18 proves the state- 
ment of the problem. 


#20: [rl := [ot s b, 0], r Z := [ct, b, e]] 

#21: r3 := [ct + R.COS(w.t), b + R.SIN(wt), 0] 


# 22 : 


"d d d 

— | rl - r2|, — | rl - r3 \ , — \r2 - r3\ 
_dt dt dt 


= [ 0 , 0 , 0 ] 


III.3 A body of mass m is moving in the horizontal plane. The initial velocity is v 0 . Find the time 
dependency of velocity and position if the force acting on the body is proportional to: 

a) its velocity, 

b) squared velocity. 


Solution: Let us write down the equations of motion: 

dv T 2 


a) 


dv T 

m — = -/v, v , 

dt 1 


b) 


m — = -k 1 v 
dt 


k x , k 2 - positive constants 


The separation of variables leads to the integral form: 



Solution of the above equations with respect to v gives the velocities. 
Let us start the Computer Algebra. 


a) We enter the first of the above equations and solve for v. To get position of the body we assign 
solution #20 to the variable v. Then we enter and simplify the equation to obtain x. 
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#24: SOLVE 


f 

V ^ 

m 

1 t 

. 

du = / 1 du, v 

kl . 

u 0 


vQ j 


kl- t/m 


#25: 

v = vQ-e 


i — i 
LX 

I 

#26: 

v := vQ-e 


t 

#27: 

X 

ii 

< 

CL 

r+ 

II 


- kl-t/m ^ 


m ■ vQ m - vO - e 

x = - 

^ kl 


kl 


b) The evaluation steps for the second force are similar. The results are displayed below: 


SOLVE 


#23 : 


k2 


1 t 

— du = J 1 du, v 

2 0 


vQ 


m- vQ 


v = 


k2-t-vQ + m j 


m- vO 


#29: v := 


k2-t-v0 + m 


#30: x = f v dt = 

0 


f m*LN(k2- t- vD + m) m-LN(m) 
x = - 


k2 


k2 


#31: v := 



f 

V 

1 




1 

t 

- kl-t/ 

#32: SOLVE 

m« 

du = 

/ 1 du, v 

= (v = vO-e 



- kl-u 

0 




vO 




t 


y 

- kl-t/m ^ 

r 


- kl-t/m 

m- vO 

m ■ vO ■ e 

#33: x = J 

vO-e dt = 

X = 


0 


kl 

kl 


Comparing the results for the considered forces we can see that in both cases the velocity tends to zero 

with time ( lim = 0 ), however the distance passed by the body when time tends to infinity only in the 

t-> o 

case of the force being proportional to velocity, for the second force (proportional to velocity squared) 
the distance tends to infinity (when evaluating the limit it is necessary to declare domains for the pa- 
rameters m , kl and k2). 

Exercise: Solve the problem for the force being the sum of both forces and discuss the asymptotics of 
solutions as dependent on values of the coefficients k\ and k 2 . 
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Titbits from Algebra and Number Theory(15) 

by Johann Wiesenbauer, Vienna 

Programming in DERIVE has always been sort of a challenge. Due to the very limited amount of pro- 
gramming constructs, it’s more like riding a wild horse rather than driving a car. At least, that’s what I 
used to say. The advent of version 5 of DfW might bring about a dramatic change on that score though 
and programming in DERIVE is sure to become a lot easier in the very near future. Hence, with one 
auspicious and one dropping eye, as it were, let’s mount that wild horse for what could be the very last 
ride. (As always the following programs refer to the latest version of DfW which is 4.1 1 at present.) 

As you may remember I wrote in my last column that STIRLINGl(n,k) is also the number of permuta- 
tions on n letters having k cycles, whereas STIRLING2(n,k) is the number of ways to partition a set 
of n elements into k disjoint subsets. Let’s call those numbers Sl(n,k) and S2(n,k), respectively. Then 
e.g. the implementation of Sl(n,k) in order to check the equality STIRLING l(n,k) = Sl(n,k) for some 
values of n and k, could be just the challenge we’re looking for. We need quite a few routines to this 
end, but as you will see, they are all interesting on their own. 

In the first place we need a routine to generate all elements of a symmetric group S n . Assuming that 
its elements are represented as the permutations of the set {l,2,...,n}, a recursive solution of this 
problem could look like this. 

SYMGROUP(n) := IF(n = 1, [[1]], SORT(APPEND(VECTOR(VECTOR( 

IN SERT ELEMENT (n, v_, n - k_ + 1), v_, SYMGROUP(n - 1)), k_, 1, n)))) 

SYMGROUP(3) - [[1, 2, 3], [1, 3, 2], [2, 1, 3], [2, 3, 1], [3, 1, 2], [3, 2, 1]] 

Next we need a routine that determines the number of cycles of a given permutation in S n . I’ll do a 
little more by providing a routine that converts a given permutation into its representation as a product 
of disjoint cycles. Actually I have already written such a program (cf. DNL #28) but the following 
version is much more straightforward and considerably shorter. 

TOCYCLES(p) := (ITERATE(IF(m_ = s_, [APPEND(a_, [c_]), [rSUBl], pSUB(r SUBl), 

DELETE ELEMENT (r_), rSUBl], [a_, APPEND(c_, [m_]), pSUBm_, SELECT(t_/= m_, t_, r_), 
sj, [a_, c_, m_, r_, sj), [a_, c_, m_, r_, sj, [[], [1], pSUBl, [2, ..., DIMENSION(p)], 1]))SUB1 

TOCYCLES([2, 3, 1, 4, 6, 5]) = [[1, 2, 3], [4], [5, 6]] 

In order to understand better how it works we replace ITERATE by ITERATES in TOCYCLES() and 
delete the SUB1 at its end thereby getting 

TOCYCLESO(p) := ITERATES(IF(m_ = s_, [APPEND(a_, [c_]), [r_SUBl], 
pSUB(r_SUBl), DELETE_ELEMENT(r_), r_SUBl], [a_, APPEND(c_, [m_]), pSUBm_, 
SELECT(t_/= m_, t_, r_), s_], [a_, c_, m_, r_, s_]), [a_, c_, m_, r_, s_], [[], [1], pSUBl, [2, 
DIMENSION(p)], 1]) 

TOCYCLESO([2, 3, 1, 4, 6, 5]) 
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[] [1] 2 [2, 3, 4,5, 6] i 

[] [i,2] 3 [3, 4,5, 6] i 

[] [1,2,3] 1 [4,5,6] 1 

[123] [4] 4 [5,6] 4 

[[1,2,3], [4]] [5] 6 [6] 5 

[[1,2,3], [4]] [5,6] 5 [] 5 

[[1,2,3], [4], [5,6]] [?] [2, 3, 1,4,6, 5] [] ? 

7 

[[1,2,3], [4], [5,6]] [?] [2, 3, 1,4,6, 5] [] ? 

7 

Now the meaning of the 5 variables a_,c_,m_,r_ ,s_ should be obvious: a_ is the current list of disjoint 
cycles. It is empty at the beginning and will contain the result we are looking for at the end of the run. 
c_ is the cycle that is currently generated. As soon as the numbers m_ and s_ coincide, where m_ is the 
image of the current end of the cycle under the permutation and s_ is its first element, this cycle is 
added to a_ and initilized with the first element of r_ which contains all elements that haven’t occurred 
so far. Eventually this r_ becomes empty thereby producing an error which is the signal for us to exit. 
(Even errors can be useful at times!) 

We have all the ingredients now to achieve our goal. 

Sl(n, k) := DIMENSION(SELECT(DIMENSION(TOCYCLES(pJ) = k, p_, SYMGROUP(n))) 

STIRLING l(n, k) := IF(k >= n, MAX(1 - k + n, 0), IF(k, 0, IF(k - 1, (n -1)!, IF(k = n - 1, 

COMB(n, 2), IF(k > n/2, (ITERATE([n_ DELETEELEMENT (v_)+ DELETEELEMENT (v_, 1 - n_ 
+ n), n_ + 1], [v_, n ], ITERATE([n_APPEND( 

DELETE ELEMENT (v_), [0]) + v_, n_ + 1], [v_, n ], ITERATE([n_APPEND(v_,[0]) + 

IN SERTELEMENT (0, v_), n_ + 1], [v_, nj, [[1], 0], n -k), 2k - n), n - k))SUBlSUBl, 
(ITERATE([n_ DELETE ELEMENT (v_, 1 - n_ + n) + DELETE ELEMENT (v_), n_ + 1], [v_, n ], 
ITERATE([n_v_ +APPEND(DELETE_ELEMENT(v_), [0]), n_+ 1], [v_, n ], 

ITERATE( [n_IN SERT ELEMENT (0, v_) + APPEND(v_, [0]), n_ + 1], [v_, n ], [[1], 0], k), n- 2-k), 
k))SUB 1 SUB 1))))) 

VECTOR(S 1 (6,k_),k_,0,6)=[0, 1 20,274,225,85, 1 5, 1 ] (72.2s) 

VECTOR(STIRLINGl(6,k_),k_,0,6)=[0, 120,274,225,85, 15,1] (0.0s) 

Although the execution times on my Pentium 166 PC speak volumes - well, it was clear from the start 
that Sl(n,k) wouldn’t be a serious competitor for STIRLING l(n,k), wasn’t it - these results show that 
our routines work. In particular, it took less than 0.015s on average to find the representation as a 
product of disjoint cycles for each permutation of S 6 , which is not so bad. 

Stirling numbers of the first and second kind appear in many areas of mathematics, in particular they 
play an important role in a number of interpolation formulae and in the calculus of finite differences. 
As I already mentioned there are also a number of applications in algebra and number theory. For 
example, you can use them to answer the question whether two given polynomial functions g(x) and 
h(x) are equal over some residue class ring mod m. After substituting f(x):=g(x)-h(x) this leads to the 
simpler question whether f(x) = 0 mod m holds for all x or not. (Note that for large n checking of all 
possible residue classes mod n might not be feasible!) Using the formulae 

: = 2jS2(n,k)x- (n>0) (*) 

k=0 


X 
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where x- : = x(x - l)(x - 2)...(x - k + 1) (k > 0) are the factorial powers and S2(n,k) are the Stir- 
ling numbers of the second kind, we can represent f(x) in the form 

f(x) = a n x- +...UJX- + a 0 

Now it can be proved that f(x) is identically zero mod m if and only if 

lc!a k = Omodm, k=0,l,...,n 
which is obviously much easier to check. 

Although the following example is a very small one (in order to save valuable space!) you will clearly 
see how this works in general. We are assuming here that the DERIVE-implementation of Stirling 
numbers of the second kind, namely 

STIRLING2(n, k) := if(n<k,0,ITERATE(APPEND(v_, [0]) + IN SERT ELEMENT (0, v_), 
v_, [1], k) VECTOR((-l) A (k - j_) j_ A n, j_, 0, k)/k!) 


has already been loaded. 

m := 2040 

5 4 3 2 

F(x) := 17 x - 170 x + 1275 x - 2890 x + 1768 x + 22440 

(a :== lim TERMS(F(x))) = [17, -170, 1275, -2890, 1768, 22440] 
x-*l 


t : 
1 
0 
0 
0 
0 
0 


= UECTOR(UECTOR(STIRLING2(n, k), k, 5, 0, -1), n, 5, 0, -1) 
10 25 15 1 0 

16 7 10 

0 13 10 

0 0 110 

0 0 0 1 0 

0 0 0 0 1 


17, 0, 680, 0, 0, 


22440] 


MOD( [17 -5? , 04?, 680 3?, 02?, 01?, 22440 0?], m) 


[ 0 , 0 , 0 , 0 , 0 , 0 ] 


Note that it may be a bit more difficult to get hold of the vector a of coefficients of f(x) in case there 
are vanishing coefficients of f(x) after the leading one. Since our construct using TERMS( ) leaves 
them out you will have to add them manually or use a different idea to compute a. 


And here another nice application of Stirling numbers before we leave this topic for good. It deals 
with computations of powers 

a 11 mod m, n=0,l,2,... (**) 


where a, k and m are fixed positive integers. Of course, we can’t beat the time DERIVE needs for the 
computation of a single power of this type, but as I will show this changes dramatically if it comes to 
the computation of a lot of such powers, say for all n< s, where s is not too small. To this end consider 
the set S of all unary vectors of a fixed length r and with nonnegative integer components < m, on 
which a unary operation f: S — > S is defined in the following way: 


D-N-L#34 


Johann Wiesenbauer: Titbits 15 


p 43 


g((x 1 ,x 2 ,...,x r )):= (XjX 2 modm,x 2 x 3 modm,...,x r _ 1 x r modm,x r ) 

i.e. x { x i+1 mod m becomes the new x i for i= 1,2,,..., r-1, whereas x r doesn’t change at all. Now let’s 
assume that we start with the vector x whose components are 

Xj := a l!S2(k,l) modm, i=0,l,2,...,k and x k+1 := 1. 

Then the first components of the sequence x,g(x), g(g(x)), g(g(g(x))), ... are exactly the numbers of 
the sequence (**) in the same order! 

Again let’s illustrate this with a small example, where a=2, k=3 and m= 1000000. The obvious way to 
compute the elements of (**) for n=0,l,2,. . .,10 would be 

[a := 2, k := 3, n := 10 6 ] 

U ECTOR (mOD ( 2 ° , m], n, 0, 10 

[1, 2, 256, 217728, 551616, 26432, 977536, 926208, 84096, 734912, 69376] 

Now let’s check our assertion by means of the following routines: 

G(x) := APPEND(VECTOR(MOD(xSUBk_xSUB(k_ + 1), m), k_, 1, k + 1), [1]) 

ITERATES(G(x), x, APPEND(VECTOR(MOD(a A (i!STIRLING2(k, i)), m), i, 0, k), [1]), 9) 

1 2 64 64 1 

2 128 4096 64 1 

256 524288 262144 64 1 

217728 953472 777216 64 1 

551616 693952 741824 64 1 

26432 248448 476736 64 1 

977536 105728 511104 64 1 

926208 3712 710656 64 1 

84096 955072 481984 64 1 

734912 422848 846976 64 1 

69376 107648 206464 64 1 

In fact, the numbers of the first column coincide with the numbers calculated above! The question 
remains whether this method is faster than the “brute force” method which computes those powers one 
by one. Well, if we compute all powers for n< s and neglect the time for the computation of the start 
vector x involving Stirling numbers of the second kind (we can certainly do this if s is large) then we 
get the matrix above by performing roughly ks multiplications. DERIVE uses internally the “square 
and multiply”-method which takes on average 1.5 k In n multiplications for one term in (**) and hence 
1.5 k In s! multiplications altogether. Using the approximation 

s 

In s! « | lntdt = s(ln s -1) 

i 

we finally get about 1.5 ks (In s -1) for the average number of multiplications, which is certainly a 
difference, if s is large! 
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The fast computation of the sequence (**) by means of Stirling numbers can be used to factor large 
integers very efficiently if certain preconditions are fulfilled and maybe I will deal with this interesting 
application another time. For now I would like to keep a promise and turn to a completely different 
topic which refers to the interesting article by Richard Schorn in DNL #33, p22-23. It goes without 
saying that I was very delighted to see a nice application of DERIVE to a problem of abstract algebra, 
which is certainly not very often the case. On the other hand, it soon became clear when skimming 
over that article that a lot of things had been done manually which should have been handed over to 
DERIVE. What follows are some ideas on that score. 


First of all, we assume that all elements of the group in question are stored in a vector 1 along with a 
symbolic name. The vector 1 typically looks like this 


e x 


1 = 


a 


x 

x-1 


b 3-x 


which means that we have chosen the names e, a,b for the functions x, x/(x-l) and 3-x, respectively. 
We now define two operations OPO(u,v) and OP(u,v) which accept only symbolic names as input, but 
the output will always be a function of x for OPO(u,v), whereas OP(u,v) additionally tries to find the 
symbolic name for the output, if one is available. 

OPO(u,v):=ITERATE((SELECT(l_SUB l=u,l_,l))SUB 1 SUB2,x,(SELECT(l_SUB l=v,l_,l))SUB 1 
SUB2,1) 

OP(u,v):=ITERATE((APPEND(SELECT(l_SUB2=o_,l_,l),[[o_]]))SUB 1 SUB 1 ,o_,OPO(u,v), 1 ) 

Note that only OPO(u,v) refers to the specific group which consists of linear rational functions in x in 
the case at issue with the composition as binary operation. At any time we can evoke the routine 
TABLE which yields the current operation table using OP(u,v). 

table:=VECTOR(VECTOR(OP(u,v),v,F SUB 1 ),u,F SUB 1 ) 

A new variable is defined by means of DEF(newvar,fl,f2) where either newvar is the composition of 
already defined variables fl and f2 (fl and f2 are supposed to be their symbolic names!) or fl is an 
expression in x defining a new linear rational function and f2 is left out. Programming buffs might 
profit a lot when looking into all the tricks I used to achieve this goal! (By the way, for the very last 
time I used here an auxiliary function like NEWL(u) to make self-assignments for 1 possible. In DfW 5 
they no longer cause problems. Good riddance!) 

NE WL(u) : =(1 : =u) 

DEF (newvar, fl ,f2):=IF(VARIABLES(fl)=[x],l SUB DIMEN SION (NE WL( APPEND(1, 

[[newvar, fl]]))), 0,1 SUB DIMENSION(NEWL(APPEND(l, [[newvar, OP(fl,f2)]])))) 

Finally I defined a routine ORD(u) that computes the order of u. 

ORDO(u):=DIMENSION(ITERATES(u,x,x))-l 
ORD(u):=ORDO((SELECT(l_SUBl=u,l_,l))SUB 1 SUB2) 

To see these routines at work, let’s set up the operation table for the same group that was considered 
by Richard Schorn in his paper. We start with the following definitions: 
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1 : = [] 

DEF(e , x) = [e,x] 

DEF | a, — — 
x-1 


a. 


x-1 


DEF(b,3-x)= [b,3-x] 

Now the vector 1 consisting of all currently defined variables has the form given above as an example. 
By evoking TABLE we get 

i a b 

x-3 


table = 


x- 2 


b 2- 


x-1 


which clearly shows how we should continue, e.g. by entering 


DEF(ab,a,b)= 


ab,- 


x-3 


x-2 


After entering DEF(ba,b,a) and DEF(aba,a,ba) in the same way (note that all these inputs must be 
simplified immediately otherwise they won’t take effect!) we finally arrive at 


table= 


e 

a 

b 

ab 

ba 

aba 

a 

e 

ab 

b 

aba 

ba 

b 

ba 

e 

aba 

a 

ab 

ab 

aba 

a 

ba 

e 

b 

ba 

b 

aba 

e 

ab 

a 

aba 

ab 

ba 

a 

b 

e 


UECTOR( ORD( 1_) , 1_, 1 1 )= [1,2, 2, 3, 3, 2] 

1 

Try it out yourself! On this occasion I want to thank Richard Schorn for many suggestions. Among 
other things he gave me the “homework” to look into the group generated by the two functions 1/x, 
and (x*(sqrt(3)+2)-l)/(x+sqrt(3)+l) which I would like to pass on to you. Have a lot of fun! 


J. W iesenbauer@tuwien.ac . at 


About Difference 1 st order Difference Equations: 

Solve the difference equations: y(x+1 ) = 1 .5 y(x) with y(0) = 5 and u(n+1 ) = 1 .5 u(n) + 2 with 
u(0) = 5: 


LIN1_DIFFERENCE 

LIN1_DIFFERENCE 



x, 0, 5 


t P 0, 5 


■ ■•(tI 


-t t + 2 
2 >3 - 4 
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Computergrafik mit DERIVE - Computer Graphics with DERIVE 

Maria Koth, University of Vienna, Austria 


3.4 Snowflake Curve and Koch-poly gons 

If you compose three Koch Curves to an equilateral triangle then you will get a starlike curve which is 
called Snowflake Curve. 

Substituting vO from the last example (Koch 
curves in a Square) by the vector ECK (n) , which 
contains the vertices of a regular polygon, the 
function KOCHPOLYGON (k, n) - similar to 
QUADRAT - produces a figure consisting of n Koch 
curves of order k over the sides of the given poly- 
gon. 

(You need K0CH(1) and the respective 
GENERATOR (x, y) from DNL#33.) 

ECK (n) : =VECTOR ( [COS (2*pi*i/n) , SIN ( 2 *pi*i/n) ] ,i,0,n) 

KOCHPOLYGON (k, n) : =VECTOR (KOCHI (k, [ (ECK (n) ) | j, 

(ECK (n) ) 4 ( j+1) ] ) , j , DIMENSION (ECK (n) ) -1) 

You can change ECK (n) to run clockwise through the polygon's vertices. 

This results in Koch curves outside of the polygon. 

ECK (n) : ^VECTOR ( [COS (2*pi*i/n ) , SIN (2*pi*i/ n) ] ,i,n,0,-l) 

It is nice to combine both kinds of curves in one figure. 




4. The Sierpinski Triangle 

The Sierpinski Triangle is also a classic example of a fractal. It was defined by the Polish mathemati- 
cian Waclaw Sierpinski in 1916. 

4.1 The construction 

You have to start with an equilateral triangle (side a) and remove its center triangle. Now you have the 
Sierpinski Triangle of order 1, which consists of three congruent equilateral triangles with sides all. 
Applying this principle of construction on each of these three triangles leads in the Sierpinski Triangle 
of order 2, etc. 

Each step of construction multiplies the number of triangles by three and bisects the length of the 
sides. So the Sierpinski Triangle of order k consists of 3 k equilateral triangles with sides a • (0.5) k . The 
Sierpinski Triangle is the set of points remaining after having applied the described function infinite 
times. 
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To obtain a DERIVE function S IERP ( k ) for creating the 
S.T of order k let us start with the triangle dO 
(side a = 1) given by the following vector: 

v0: = [[0,0],[l,0], [l/2,SQRT(3)/2], [0,0]] 

The Sierpinski Triangle(l) consists of this initial triangle 
vO together with triangle vl, which connects the mid- 
points of vO's sides. 



v3 



vl : = [ 0 . 5* (vO 1 1+vO 1 2 ) ,0.5* (v0|2+v0|3) ,0.5* (vO j+vO |3) ,0.5* (vO jl+v0 |2) ] 


Proceeding to the Sierpinski Triangle(2) we have to add three new triangles, which are obtained apply- 
ing three centric stretchings onto triangle dl with vl, v2 and v3 as streching centers and 0.5 as stretch- 
ing factor. 

Applying these stretchings onto each of the newly created triangles we can find 9 triangles which form 
together with the 1 + 1 + 3 = 5 existing triangles the Sierpinski Triangle of order 3. In the same way 
we apply the stretchings onto the last created 3 k l triangles when we proceed from Sierpinski Trian- 
gle^) to Sierpinski Triangle(k+1). 

STRECK ( i , v) :=VECTOR( VECTOR (0.5* (vOji+vjj jk) , k, 1,4) , j,l,DIM(v) ) 

Being v the vector of the triangles created in the last step, the STRECK (i, v) -statement causes the 
stretching of all these triangles by the factor 0.5 with the i-th vertex of vO be done. SIERP (n) pro- 
duces a Sierpinski Triangle of order n. 

SIERP (n) : =IF (n=l , vO , APPEND ( [vO] , APPEND ( ITERATES (APPEND ( STRECK ( 1 , v) , 
STRECK (2, v) , STRECK ( 3, v) ) ,v, [vl] ,n-l) ) ) ) 



SIERP (8) 
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4.2 Sierpinski Polygons 

Construction of a Sierpinski Triangle can also be explained as follows: 

An equilateral triangle is substituted by three congruent equilateral triangles which 

• are parallel to the sides of the initial triangle, 

• have one vertex with the initial triangle in common, 

• have side lengths such that two neighbouring triangles have one point in common. 
This principle of construction can easily be transfered on regular n-gons with n > 5. 




The small n-gons of the next generation are 
created by an appropriate stretching with one 
vertex of the parent-gon as center. The 
stretching factor f depends on the kind of 
polygon. If it is a triangle, then we need 
three stretchings with factor f = 0.5 - as we 
have seen -, if it is a regular hexagon then 
we need six stretchings with factor f = 1/3 
etc. 


If we start with a regular n-gon with side length 1 = 2/ + x, then we will call the side length of one 
’’child”/ The sketch above shows how to find /in an easy way: 


2 1 + cos^ 
v n 

For the DERIVE procedure we first define / and the coordinates vO of the initial n-gon: 

[f :=!/ (2+2*COS (2*pi/n) ) , vO :=VECTOR ( [COS (i) , SIN (i) ] , i , 0 , 2 *pi , 2 *pi/n) ] 

Like working with the Sierpinski Triangle the function TEIL ( v, i ) generates the i-th partial polygon 
of the parent-polygon v. 

TEIL (v, i) :=VECTOR ( ( 1— f ) *v SUB i+f*v SUB j,j,l,n+l) 

TEILALL(v) creates the complete next generation of figure v 
TEILALL (v) : = VECTOR (TEIL (v, i) , i, n) 

As the Sierpinski polygons of order k - v = SIERPP(k) - consists of n k polygons v i? we have to apply 
TEILALL(vi) onto each of these polygons Vi. The recursive statement SIERPP(k) creates the polygon 
of order k - after having defined n! ! 

NECHST(v) : =IF (DIMENSION ( v) =1 , TEILALL ( v SUB 
1) r APPEND (NECHST ( DELETE_ELEMENT 

(v, DIMENSION (v) ) ) , TEILALL (v SUB DIMENSION (v) )) ) 


SIERPP (k) :=IF (k=0, [vO] , NECHST (SIERPP (k-1) ) ) 
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5. The Chaos Game 

1985 the English mathematician Michael Barnsley presented a new interesting method to create fractal 
pictures. He called this method “Chaos Game”. 


5.1 The Chaos Game for the Sierpmski Triangle 



Start with an equilateral triangle with vertices 
vl, v2 and v3. Then pick any arbitrary point zO 
anywhere in the plane. Now, randomly select 
one of the three vertices, say v2. The midpoint of 
the segment v2,z0 gives the next point zl. Select 
- again randomly - one of the vertices, say vl — » 
midpoint of segment (vl,zl) gives z3, etc. 

In DERIVE we can describe this algorithm using 
the RANDOM(n)-function, which returns for 
each n e N and n > 0 a randomly chosen element 
of {0,1,2, ,n-l}. 


So l+RANDOM(3) gives a randomly chosen element of {1,2,3}. The rule to obtain the sequence of 
points is given by: 

Z n +1 = 0.5(z n + Vi +RANDOM(3))- 

The DERIVE code consists of only two lines. We define the initial triangle v together with the starting 
point zO, e.g. [1,1]. CHAOS(n) produces n+1 points of the Chaos Game: 

[v:=[[0,0],[l,0], [0.5,0.5*SQRT(3) ] ] , zO := [1, 1] ] 
random ( 0 ) = 471223100 


CHAOS (n) : =ITERATES (0.5* (z+vj (l+RANDOM(3) ) ) , z, z0,n) 
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It is important to plot the points in the "Discrete” Mode. What might be the result for large numbers n? 
We would expect a random distribution of points within the initial triangle (- trivial for a starting point 
within the circumference, and easy to explain that if starting from outside we will come closer to the 
triangle and at a certain moment a point will be captured inside and all the following points will stay 
inside). But the distribution shows - although unpredictible for each point - a predictible pattern: 

CHAOS (10) 

CHAOS (100) 

CHAOS (1000) 

CHAOS (10000) 

Very clear appears the Sierpinski Tri- 
angle out from the shadow of points! ! 



(You can find a TI-92 program for the Chaos Game in "Mastering the TI-92" from Rich, Rose and 
Gilligan, GILMAR 1996 - DNL#22, Josef.) 

The principle of construction contains the same stretching as the construction of the Sierpinski Trian- 
gle in 4.1. If you start the Chaos Game with a point zO which is an element of the Sierpinski Triangle 
the sequence of the following points will all be points of the Sierpinski Triangle If you choose another 
point then the points of the sequence approximate the Sierpinski Triangle Let x be a point of the Sier- 
pinski Triangle and 5 > 0, then there exists an element of the sequence which is lying in the 5- 
surrounding of x. [Peitgen]. 


5.2 More Chaos Games 

In 4.2 we transferee! the Sierpinski Triangle Construction onto polygons. We can do the same with the 
Chaos Game. 

[v6:=VECTOR ( [COS (i) , SIN (i) ] , i, 0, 5*pi/3+0 . l,pi/3) , f :=l/3, zO := [1, 1] ] 
CHAOS 6 (n) : = ITERATES (f*z+ (1-f) *v6j. (l+RANDOM(6) ) , z, z0,n) 

CHAOS 6 (6000) 

Needs only 3 seconds for being plotted (DERIVE 6). 

In 1999 46.5 seconds were necessary! 

(I added the 0.1 within the VECTOR command 
that “Approximate before plotting” can be used in 
the 2D-Plot Window, Josef) 
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Working with Maria Koth’s file from 1999 I had the idea to try a more general function includ- 
ing the kind of polygon and the number of iterations as function parameters: chaosp(k,n) was 
the result: 


f_(k) := 


f f 2, "Y 

2-1+ COS 

l l k 1 


f 2 (k - 1) 2 ■ it 

v_(k) := VECTOR [CQS(i), 5IN(i)] P i, 0, + 0.1, 


k k k j 

chaosp(k, n) := ITERATES(f_(k) ■ z + (1 - f_(k)) ■ (v_(k)) , z 9 z0 ? n) 

1 + RANDOM(k) 


chaosp(S, 4000) 



chaospClZ, 4000) 
chaospClZ, 4000) 
cfiaosp(12, 4000) 

Three plots superimposed (12 000 points) 


References: 

Endl, K. Kreative Computergrafik. VDI Verlag, Diisseldorf 1986. 

Hermann, D. Algorithmen fur Chaos und Fraktale. Addison- Wesley, Bonn 1994. 
Koth, M. Einfuhrung in das Arbeiten mit DERIVE. Pl-Wien, Wien 1994 
Peitgen, H.O. u.a. Fractals for the Classroom , Part 1. Springer, New York 1992. 
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It is not difficult to present the Chaos Game with 
TI-NspireCAS. 

I thought that I would need a LUA program in 
order to produce a plot displaying 2000 points or 
more. Then I was surprised that - at least on the 
PC - lists of so many elements can be used for 
producing a scatter plot. 

The parameters for the Game are entered in a 
Notes page. 


chaos 

Define cliaos()= 
Prgm 

Local v t nst 

x/: = [sM f l] } 

yl={st[ 1.2]} 

nst. =st 

For i, 1 ,n 

v. =randlnt(l j 3 ) 

n$£+£r\v 1 
_ L J 


2 

x/:=apprax(augment(x/, { 
yl. =approx (augment^/, { 
End For 
End Prgm 

^[u]})j 

^[ 1 , 2 ]})) 


CHAOS GAME 

The vertices of the given equilateral triangle: 



-9 9 


-3^3 -3 

ti : = 

— ,-3— ,-3:0,6 
J3 J3 


3' J3 -3 
0 6 


Enter a starting point st with -10 £ x £ 10 and -6 £ y £ 6: st:=[2,2] ► [2 2 ] 
Enter a random seed number: RandSeed 2111 ► Done 
Enter the number of iterations: n:=2000 ► 2000 
Run chaos(): chaosQ * Done 

Lists xl and yl must be displayed in a Scatter Plot). 
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CHAOS POLYGONS 

Enter k for a k-polygon: k:=8 ► 8 

Enter a starting point st with -10 £ x £ 10 and -6 £ y £ 6: st:=[l,l] 
Enter a random seed number: RandSeed 1111 
Enter the number of iterations: |n:=3^0d 
Run chaosp(): chaosp() 

Lists xl and yl must be displayed in a Scatter Plot. 


The Chaospolygons need one function and one 
program. The chaotic octagon is plotted with two 
different attributes. 



chaosp 

Define chaosp ()= 
Prgm 

Local v f nst 
nst. =st 



For i, 1 ,n 
■: =randlnt( l,fr) 

nst =appnox \f- nst+ ( 1 ■ p ( k) [■ v] ) 
xl. =augment(x/, { nst _ 1, 1 
: =au gment(v'4 { nst 1, 2 
End For 
End Prgm 


11 


P 


Define p(/r)= 
Func 
Local m f i 
m\ =newMat(/c J 2) 
For i , l } k 


-W:- 


2 ' 71 f \\ J 2 JT f \ 

cos[ ■ (z-lM sin ■ (j-l) 

k I \ k 


End For 
nt\ 

End Func 
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An Investigation in Calculus using the TI-92 

by Michael Meagher, Vienna International School 


We start with the problem as follows. Given 
the graphs of y = x, y =1 and x=l: 


rriTTOT Y F3 Y FH Y FEt Y FfiT YF? a 
|ZQQn|Trace|ReGraph|Maih|Drau|T y 


We can therefore see that in general the ratio 
of Area A : Area B is n : 1. 

We would like to develop the problem in a 
number of ways, in particular by changing the 
limits of integration. We would like to use the 
TI-92 for this problem and so, as a first step, 
we will check our work thus far on the TI-92 
by defining some functions to calculate the 
areas A and B. 


We then ask the simple question: What is the 
ratio of the Area A to the Area B? By straight- 
forward calculation of the areas of the two 
triangles we get the answer 1:1. 

The problem starts to become more interesting 
if we consider the graphs y=x A 2, y=l and x=l. 


| Define below<n>=.f <x A n J .x J .0 J .l>| 


Whenever a function is defined on the TI-92, it 
is always a good idea to test the function on a 
few known results: 


■ below(2) 



1^3 

■ below(3) 



1^4 


MAIN F: 

AD EXACT 

FUNC iB. 



All seems well. We define above(n) and test 
some results: 



The area B is calculated, by hand, as x 

JO 


which is 1 / 3 and, since the area of A is calcu- 
lated by subtracting the area of B from a unit 
square, we get that the area of A is 2/3. The 
ratio of the areas is therefore 2:1. 


And now for the interesting result: 


abo^e(n> 

belou(n) 


Following the same steps for y = x A 3 we find 
that the area is 3 : 1. We conjecture at this point 
that there is a relationship between the index 
of x in the curve and the ratio of Area A : Area 
B. This ratio can be calculated as follows: 


What???? 

Look at all the examples we have seen so far: 
n= 2, 3, 4, 50. All of these are positive integers. 
Of course the original question was posed 
with n implicitly meant to be a positive inte- 
ger. A moment's thought about the graph of y 
= x A n with x a negative integer reveals that 
there is no area enclosed under the graph be- 
tween 0 and 1 as the following graph of x A (-2) 
suggests: 


n + 1 


n + 1 










This is a good example of how the TI-92 pro- 
motes exact mathematical thinking. Our casual 
definition of the function below (n) was not 
quite what we meant but the TI-92 forces us to 
think more clearly. 



We need to modify our definition of below(n) In general the value of f(x) = x A n at x = 2 is 

using the "with" ( | ) operator: 2 A n, so the rectangle we need to consider for 

the total area is 2 * 2 A n or 2 A (n+l). So procee- 



Now we go back to our question: 


Define below(n) = g(x n Jdx|n>0 


1 Define aboue(n) =2 n + * - below(n) 


Lne above <n>=2 A <n+l>-belou<n> 





A reasonable first extension to this problem is 
to ask what happens if we integrate from 0 to 
2. Before we start again we want to delete the 
functions below and above: 



We start again from the basics and can gener- 
ate the first results by hand or with the TI-92 
as follows: 


Extension 2 

So if this property holds while integrating 
from 0 to 1 and integrating from 0 to 2, will it 
work from 0 to b, an arbitrary parameter? 



belou<n J .b>=J<x A n J .x J .0 J .b> I n>Ol 


We check this definition for known values: 



Not very promising, but of course the total of 
area A and area B is no longer a square. The 
value of the function f(x) = x A 2 at x = 2 is 4 and 
so we are looking at a two by four rectangle 
with area 8. 


The value of the function f(x) = x A n at x = b is 
b A n. The rectangle under consideration is then 
b*b A n or b A (n+l) : 
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I 


J 


■ Define aboveCn , b) = b n + * - belouCn , b) 

Done 


... bove<n J .b>=b A <n+l>-belQw<n 


IMflIH m EHfltT FUH'. 1WW 


And we look for the general result: 


abowe(n ? b) 

below(n ? b) ^ 


above <n J .b>/belou<n J .b> 


IMflIH FirtD EKflCT FUNC 30/30 I 

And so the property holds on the integral from 
0 to an arbitrary parameter b. 


Extension 3 

The next question that arises is: if it works for 

0 to 1 and it works for 0 to 2, does it work from 

1 to 2? 



I MAIN RAD EiiflCT FUNC 

As before we delete below and above: 

n -r 1 

I" DelUar abowe, below Done| 


delvar above ^belowl 


IMflIH m EHfltT FUliL 3WW 

To find the new area B we take: 


■ 1 1 (x 2 )dx ■+ below 7x3 


j <x A 2 ..xA ,2>-frbelow| 


IMflIH m EKftCT FUHt I 

As before the total area is a 2 by 4 rectangle 
with area 8 and there is a unit square left un- 
shaded. So the area Ais8-l-7/3 


14x3 

2 


above/beloul 


IMflIH m EKfltT FUHt I 

And again the property seems to hold in this 
case. 

We now turn our mind to the most general 
result which is for the curve x A n, considered 
between arbitrary parameters a and b with a < 
b: 


■S - 1 - below + abowe 
abowe 
below 


The values of the function f(x) = x A n at a and 
at b are a A n and b A n respectively. The total 
area of the shape is b*b A n or b A (n+l) and the 
unshaded part is a*a A n or a A (n+l). Therefore: 


■ Define abovetn , a , b) =b n + 1 - a n + 1 - be ► 

Done 


... b A <n+l>-a A <n+l>-below<n ..a^b^ 


IMflIH m EHfl’.T FUliL 

Once more we check the definitions for some 
known values: 


■ belou(2 ? 1 „ 2) 

■ abooe(2 , 1 , 2) 
abooe(2 ? 1 ? 2) 
below(2 j 1 , 2) 


above <2 ,.1 J .2>/below<2 ,.1 ,.2> 


IMflIH FirtD EKflCT FUNC 30/30 


7x3 

14x3 

2 


Finally, the general result: 


a abooe(n ? a ? b) 
below(n j a ? b) 


above<n j.aj.b^/belowCn J .a J .b> 


IMflIH Kfln FHAfT FIIMf Wifi 

Extension 4 

This entire investigation can be expanded to 
the question of volumes of revolution. It works 
but doesn't have the same answer as with the 
area and is dependent on the axis of revolu- 
tion. 


Extension 5 

To return finally to the question of a negative 
value of n. The property does, in fact hold for 
values of n so long as the limits of integration 
are either both positive or both negative. In 
other words the problem is if one of the end- 
points of the integral is 0 or the integration is 
across 0. 

To my surprise however, when we define the 
functions below(n,a,b) and above(n,a,b) with 
no restriction on n we get 



FIt 

Algebra 

fFV^-1 

Calc 

" FHt " 

Other 

" FE 

PrgmlO 

f Ffi ) 

Clear a-z... 

■Define belowCn n a n b) = 

^ [x n Jdx Done 


■ Define above(n , a , b) =b n + 1 - a n + 1 - be ► 

Done 


a abooe^n p a ? b) 
below^n , a ? b) 
u abooe^ ~2 ? a ? b) 
beloiX ~2 p a p b) 


[above <~2 J .a J .b>/<belou<~2 ,a,b>> 


MAIN RAD EHACT FUNC H/30 
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As I mentioned on page 24 I wanted to solve linear difference equations of 1 st and 2 nd order 
with the TI-NspireCAS. For this purpose I inspected the respective DERIVE functions in the 
utility file provided by DERIVE. 


LIN2_CCF_H0RLAUX(b , x) := 

If b < 0 

#1: - COS(jt-x) ■ (-b) A x 

b A x 

LIN2_CCF_HQM(p , q, d r x) := 

If d = 0 

[(- p/2) A x, x-(- p/2) A x] 

#2: If d > 0 

[LIN2_CCF_HOMJ5.UX((-p Wd)/2 r x) . LIN2_CCF_HDMJVUIX((-p - -/d)/2, x)]| 

C(p A 2 - d)/4) A (x/2) ■ [SIN(x-ATAN(V(-d) , -p)) , C0SCx.ATAN(7C-d) , -p))] 

(( P A 2 - d)/4) A (x/2).[SINtx.ATAN(/C-d), -p)) p GDS(x-ATAHCJ(-cO , -p))J 

[- ELEMENTS, 2), ELEMENT(h , 1)] 

LIN2_CCF_PARTIC(h , r, x) := h-E r- Tim 

#3: x x-*x + 1 DET(Th P Tim hi) 

L x-»x + 1 J 

#4: LIN2_CCF_AUX(h r r P x r c) := c-h + LIN2_CCF_PARTIC(h r r P x) 

2 

#5 : LIN2_CCF(p P q, r P x P cl P c2) := LIN2_CCF_A.UX(LIN2_CCF_H0M(p P q f p - 4-q P x) p r P x P [cl, c2]} 


Solving one example (stepwise): y(x+2) + 6 y(x+1) + 9 y(x) = 10 

LIN2_CCF_HQM(6, 9, 0, x) = [(-3) , x-(-3) ] 


2-x + 1 2 2 -x + 1 

DET([h P lim hi) = 3 ■(! - 2-CQS(rr-x) } - 2-3 ■ L -SIN&r-x) -CQS(rr -x) 

L x-*x + 1 J 

[- ELEMENTfh P 2} p ELEMENT(h P 1)] F -x - 1 -x -x - 1 -xl 

= 1,3 -x-(-l) p - 3 ■(-!) J 

DETtTh p lim hi) 

L x-*x + 1 J 

[- ELEMENT [h P 2) , ELEMENT(h P 1) ] [ -x - 2 -x -x - 2 -xl 

10- lim =1-10-3 -Cx + l)-(-l) p 10-3 -(-1) J 

x-*x + 1 DET CIT h p lim hi) 

L x-*x -f 1 J 



[- ELEMENT (h p 2) p ELEMENT Ch p 1)] 

1 10- lim 

x x-*x -if 1 DET(Fh P lim hi) 

L x-»x + 1 J 

[- ELEMENT(h P 2), ELEMENT (h p 1) ] 

h-Z 10- lim 

x x-»x + 1 DET C|~ h p lim hj) 


-x - 1 -x 

5-3 ■ (4-x + 3) -(-1) 


S 


5 


S 


-x - 
5-3 


1 -x 

■C-1) 


2 


LIN2_CCF(6 P 9 P 10 p x P cl P c2) 


5 x 

+ (c2-x + cl)- (-3) 

a 


5 x 

I check the result: y(x) — + (c2-x + cl). (-3) 

S 


y(x + 2) + 6 ■ y(x + 1) + 9-y(x) = 10 
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Function l2chom(p,q,x) solves the homogeneous difference equation, which is no problem. 
It is the equivalent to DERIVE’s LIN2_CCF_HOM. I had to define a special ATAN-function 
with two arguments: “ATAN(y, x) is the angle of the point (x, y) in the x-y plane measured 
counterclockwise from the positive x-axis", (DERIVE Online Help). 


& 12chom( 6, 9,0.x) [(l) x - 3 X (-l) x x-3 x ] 

I2cp 1/5 

Define 12cp(X?’,x)= 

Func 

Local dt,dtl 

dt =det J | au gment/ft T , | lim (&) 1 T U T 

A:=[(-l) X -3 x (-l) x - x- 3 X ] 

[cos(jtx)- 3 x +sin(jr x)- 3 X i x cos(jtx)- 3 x +x- sin(jrx)- 3 X i] 

© the determinant dt: 

11 1 \x->x+l )]), 

m - ii m (’-W 1 - 2 ! “t 1 - 1 !! 

x^x + l l * ! 

!2 c p{h, 10,x) ( 3 _ 6 . ( cos ( n . x ))2). 9 x_ 6 . sin{jt . x ) CQS ( n . x ) j*. t 

© the lim dtl multiplied by r (= lo): 

adiff{dtl,x) 

©dotP (a, adiffidt l } x ) ) 

End Func 

► 

I2cp(h, 10,x) 

-10- (x+l)- cos (it x) 3 x 10- (x+l)- sin(jr x) 3 x m 10- cos(jr x) 3 x 

|_ 9 9 9 

©( the indefinite sum = the antidifference of the above expression): 


I2cp{h, 10,x) 
1 

sc-1 

-lo - 2>i>- cos (^' n )' ^ 

k= 1 -5- cosier x) 3 x 5 

u 

0! 

> 

1 CO 

1 1_— ^ 

O’. 


The hard problem was - and still is - to im- 
plement the indefinite sum (the antidiffer- 
ence) for not only polynomial functions. 

Albert Rich admits in the DERIVE Online 
Help that it is not easy to find the antidiffer- 
ence in many cases: 

“Even when such an antidifference exists, 
there is no known method that is guaranteed 
to find it.” 

Below is a 1 st order difference equation: 

y(m+1) = m-y(m) + (m+1)!; y(1) = 2 solved 
with DERIVE and TI-NspireCAS: 


adij^P' + l,x) 

x 2 + 7 x ^ 
3 2 6 

T, 9 

t v x x 7x 

m= + — i 

3 2 6 

Per tig 

i x +l)-Ax) 

x 2 + l 

adij[ 10- 3" x “ 2 - (-l) x ,x) 

-S- cos (tc ■ x) 3 x 5 
6 IS 

j \ -5- cqs(jt x) 3 x 5 

Kx): 

6 IS 

Peril g 

/(x+l)-/(x) 

10- cos(ji- x)- 3 x 
9 


LINl_DIFFERENCE(m, (m + 1)!, m, 1, 2) 


2 


m 


2 



+ 


1 


■Cm - I)' 


lifi {m + 1) ! , m, 1, 2 ) 


(ffz^+ffz+2) 


m\ 


2- m 


See two more 1 st order difference equations on page 45. 


